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ABSTRACT 



Moored current meter observations were made along the 
central California continental slope from Pt. Piedras Blancas 
to the Farallon Islands. The study area covered an alongshore 
distance of approximately 290 km and examined the data from 
off Pt . Piedras Blancas (P4) , Pt . Sur (P2 and P3) , Monterey 
Bay (MBl and MB2) and the Farallon Islands (D and E) . Time 
and frequency domain analyses were performed on three time 
segments that included data from combinations of the above 
mooring locations based on common time periods and depths (350 
m - 500 m) . Segment 1 revealed a mean poleward flow 
attributed to the California Undercurrent at all moorings. 
Segment 2 had a similar mean poleward flow but also exhibited 
an equatorward reversal at periods between 19.5 and 58.5 days 
at P2 and MB2 that was less apparent at P3 further offshore 
and appeared to be coastally trapped. Observed wavelengths 
were compared with simple wave theory indicating that the 
observed signal may have resulted from a combination of two 
theories i.e. a coastally trapped wave. The complex bottom 
topography between Pt . Sur and the Monterey Bay prohibited 
more rigorous comparison between theory and the observations. 
Satellite sea surface temperature imagery during Segment 3 
showed an anti -cyclonic meander with a 65 km radius outside of 
the Monterey Bay. Current and temperature records at P2 and 
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MB2 indicated that the surface feature's position varied and 
influenced the currents at depths down to 500 m. Time domain 
empirical orthogonal functions were calculated for alongshore 
and across -shore components separately. The first two 
empirical modes accounted for between 81% and 86% of the 

alongshore variance in all segments and were attributed to the 
mean California Undercurrent influence (mode 1) and to the 

deviations from the mean state (mode 2) . The first two 

across- shore modes explained between 66% and 81% of the 

variance and consistently demonstrated opposing current 
structures between the Pt. Sur and Monterey Bay moorings, 
likely due to the Monterey Bay meander influence. 
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I . INTRODUCTION 



A. EXPERIMENT BACKGROUND 

The paper's objective is to examine and describe the 
alongshore variability of sub- thermocline currents (350 m - 
500 m) over the California continental slope. The 
investigation employs conventional frequency domain (spectral) 
analysis as well as time domain empirical orthogonal function 
(EOF) analysis to provide a detailed picture of three 
different time segments. The analysis incorporated current 
meter data from several moorings between Point Sur and the 
Farallon Islands along the California coast (Figure 1) . The 
study encompassed an alongshore distance of approximately 290 
km. All moorings were deployed on the middle slope between 
the 800 m and 2000 m isobaths. Mooring P4, the furthest south 
of all moorings was located off Point Piedras Blancas, 
approximately 100 km to the south of Point Sur. Two moorings, 
P2 and P3 were positioned off Point Sur. Mooring P3 , provided 
data from the lower continental slope (1800 m) while P2, 
approximately 26 km east of P3 , was positioned over the upper 
slope (800 m) . MB2, located 41 km to the north of P2, was 
positioned outside of the Monterey Bay (1979 m) while MBl, 
approximately 31 km to the east of MB2 was inside the Monterey 
Bay over the axis of the Monterey Canyon. Positioned 
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Figure 1. Current Meter Mooring Locations! A map of the 
central California coastal region. Moorings are 
shown as solid circles with alphanumeric labels. 

approximately 111 km to the north of MB2 , mooring D marked the 
southern boundary of the Farallon Islands data set. Finally, 
mooring E was located 32 km north of D in the vicinity of the 
Farallon Islands outside San Francisco Bay. (Figure 1) 

The analysis of alongshore propagation required common 
time segments from each location. Various organizations 
collected data for different purposes, limiting the number of 
common time segments available. Instruments at depths between 
350 m and 500 m provided the most continuous alongshore data 
records so the analysis focused on these depths. The 
experiment used three segments of 110, 199 and 177 days that 
spanned an overall time period from December 1989 to November 
1991 (Figure 2) . 
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Figure 2. Current Meter Time Line: Common time periods are 

labeled as Segments 1, 2 and 3 and are shown by 
different cross-hatching styles. 

The remainder of this chapter provides background 
information on the California Current System and describes 
various deviations from the mean flow observed there. A 
discussion on coastal trapped waves is presented here to 
include the derivation of the dispersion relationship for 
internal Kelvin waves. The second and third chapters review 
data collection and analysis techniques respectively. The 
fourth chapter presents the results and is followed by a 
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discussion chapter that applies theory to the observations. 
The final chapter gives the conclusions with recommendations 
for future research. 

B. THE CALIFORNIA CURRENT SYSTEM 
1. The Mean State 

The data were collected from within the California 
Current System (CCS) , the eastern boundary current of the 
North Pacific gyre. This sub- tropical gyre consists of an 
anti -cyclonic surface flow bounded by the westward flowing 
North Equatorial Current (20° N) to the south, the poleward 
flowing Kuroshio current to the west and the West Wind Drift 
to the north. When the West Wind Drift strikes the North 
American continent, it splits to the north to feed the 
poleward flowing Alaskan current and to the south to form the 
California Current. (Wooster and Reid, 1963) 

Large scale atmospheric forcing creates this anti- 
cyclonic circulation. During the Coastal Ocean Dynamics 
Experiments (CODE 1 in 1981 and CODE 2 in 1982) , Halliwell and 
Allen (1987) attribute persistent , summertime equatorward wind 
stress to the interaction of the North Pacific subtropical 
high with the persistent low over the southwestern United 
States. The pressure differences between these two stationary 
systems are largest during the sxammer and hence equatorward 
wind stress is also largest in summer (Huyer, 1983) . Huyer 
(1983) describes the southwestward migration of the 
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subtropical high to about 28° N in February. Halliwell and 
Allen (1987) attribute winter time poleward wind stress to 
northward propagating cyclones and anti - cyclones . In general, 
atmospheric forcing drives an equatorward surface flow during 
summer which is often described as slow, broad and shallow 
(Wooster and Reid 1963) . When these equatorward winds relax, 
the current reverses within 150 km of the coast (Winant et 
al . , 1987; Chelton et al . , 1988) . This winter reversal is 
known as the California Countercurrent or the Davidson Current 
north of Point Conception (Chelton, 1984) . 

Since this paper focused on the 350 m to 500 m depth 
interval, the poleward flowing California Undercurrent (CUC) 
influenced the data set more than any other feature of the 
CCS. Hickey (1979) observed the CUC at a depth of 300 m off 
southern Washington in 1972. Chelton (1984) also observed 
this flow off central California and found an annual variation 
at Point Sur, in contrast with the semi-annual variation in 
the northern CCS such as that observed by Hickey (1979) . 
Wickham et al . (1987) observed a "jet-like core" with velocity 
greater than 15 cm/s confined to within 30 km of the coast 
near Point Sur. Like Chelton (1984) , they also found that the 
CUC varied annually with maximum poleward currents during May 
and June of 1979 and 1980. Chelton et al . (1988) found 

poleward flowing currents from the surface down to 400 m as 
far as 100 km beyond the central California shelf break during 
the summer of 1981 and 1984. They attributed this unusual 
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occurrence to an extended period of anomalously weak 
summertime equatorward winds. Such poleward flow had been 
previously observed only during winter as the Davidson 
Current . 

Hickey (1989) studied the shelf and slope undercurrent 
off Washington and the Southern California Bight, nearly 
encompassing the entire CCS. Of particular relevance to this 
paper, Hickey discussed interruptions or anomalous cases where 
the generally poleward flowing undercurrent became 
equatorward. For example, in winter off Washington, 
alongshore wind stress and pressure gradient were opposite in 
direction and reversed from the summer case creating a 
wintertime equatorward undercurrent off the Washington slope. 
Additionally, Huyer et al . , (1984) related equatorward flow 
over the slope to the presence of an eddy. Equatorward 
deviations likewise occurred in this data set and will be 
discussed on a case by case basis. 

2 . The Barocllnlc Equatorward Jet 

Contrary to Wooster and Reid's (1963) California 
surface current description, data from the Coastal Transition 
Zone (CTZ) experiment revealed a strong baroclinic jet (core 
speed > 50 cm s'^) between 37° N and 38° N (Ramp et al., 1991; 
Huyer et al., 1991; Kosro et al . , 1991). Huyer et al . (1991) 
superimposed acoustic doppler current profiler (ADCP) vectors 
onto satellite imagery and clearly showed a highly energetic 
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meandering surface flow. This jet displaced isopycnals to 
depths greater than 500 m (Ramp et al . , 1991) . Daggett (1989) 
combined CTD and ADCP data to describe an equatorward 
meandering jet off Point Arena, California during CTZ. 
Equatorward velocity in excess of 25 cm sec'^ was observed in 
a shallow (50 m) core with weaker equatorward flow extending 
down to 500 m. Additionally, in May 1988 Tisch (1990) 
observed an equatorward jet with velocity exceeding 25 cm 
sec'^ off Point Sur down to 200 m. This flow was attributed 
to strong (-2.0 dyne cm"^) wind stress present during the week 
prior to hydrographic observations. In November 1988, 
equatorward geostrophic velocity was observed from the surface 
to depths exceeding 600 m (Tisch 1990) . This paper provides 
additional evidence that an onshore meandering of this 
equatorward jet can interrupt the otherwise persistent 
poleward flow at depths down to 500 m (Chapter IV) . 

3 . Low Frequency Current Fluctuations 

In addition to the baroclinic jet interaction, several 
researchers have observed more periodic current fluctuations 
at periods from 5 to 40 days. Kundu, Allen and Smith (1975) 
compared dynamic and empirical modes of the velocity field 
near the Oregon coast. Although their observations were 
limited to the shelf (100 m) , they found high mutual coherence 
between the wind stress and barotropic u and v modes at 
periods between 7 and 16 days. Battisti and Hickey (1984) 
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found significant coherence between two synthetic sub- surface 
pressure (SSP) levels over the entire low frequency spectrum 
band for coastally trapped waves (CTW) with periods greater 
than 5 days. They calculated the of f shore-vertical structure 
and dispersion relationship of a free CTW and compared model 
results with observations. Again, these results were largely 
confined to the surface flow. 

Romea and Allen (1983) used an f -plane model to study 
vertically propagating, coastal trapped, internal Kelvin waves 
(IKW) . They paid particular attention to low frequency 
behavior with periods between 20 and 40 days, using three 
different types of wind forcing: an impulse wind at a point, 
a steady wind at a point, and a travelling wind. Their model 
showed that poleward travelling winds forced vertically 
propagating IKW while equatorward travelling winds forced 
waves that became surface trapped. They addressed the 
circumstances under which waves may propagate vertically 
without dissipation and their results will be compared to the 
observations in this paper in greater detail. Vertical IKW 
propagation must be considered here as it provides a viable 
means by which a remotely forced wave could propagate between 
Point Sur and the Monterey Bay. Denbo and Allen (1987) , 
Chapman (1987) and Lopez and Clarke (1989) found that remote 
atmospheric forcing to the south was an important source of 
current variability on the western U.S. continental shelf. 
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C. COASTAL TRAPPED WAVE THEORY 

A coherent, low frequency fluctuation (T = 19-29 days) 
persisted during the second time segment of this paper's data 
set. The dispersion relationships for simple internal Kelvin 
waves and barotropic continental shelf waves are presented as 
possible points of comparison for these observations. 

By using the inviscid momentum equations with the 
continuity equation, a series of wave- like solutions whose 
amplitudes decay exponentially with distance from a straight 
coastline can be produced. These solutions are generally 
referred to as "coastal trapped waves" (CTW) (Buyer, 1990) . 
Buyer (1990) described two extreme types of CTW. On one 
extreme is the continental shelf wave (CSW) which assumes a 
homogeneous fluid and relies solely on a sloping bottom for 
trapping. This wave is also referred to as a topographic 
Rossby wave and has the following dispersion relationship 
(Gill 1982) : 

0 , = 

(*2 . » A . 2 ) 

where co is the wave frequency in radians sec'^, f is the 
Coriolis parameter, k and 1 are the across -shore and 
alongshore wave numbers respectively and X is the bottom 

1 hh 

slope parameter = where x is the distance offshore and 

B is the bottom depth. The derivation of this relationship 
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assumes a depth which increases exponentially with distance 
from the coast. 

On the other extreme, Kelvin wave theory assumes a flat 
bottom with a vertical boundary to trap waves of sufficiently 
low frequencies i.e., a)<< f. Internal Kelvin waves (IKW) 

represent the higher order modes of the Kelvin wave solution 
and depend only on the fluid's density stratification (Huyer, 
1990) . The IKW dispersion relationship is derived below 
following Gill (1982): 

Applying a vertical boundary at x = 0, all u (across - 
shore) components equal zero (i.e., current may not flow into 
the wall) and the following governing equations result: 

-fv = -A 

p ax 

dv _ __1 dp 
dC p dy 



( 1 ) 



( 2 ) 



dw _ __i_^ _ p'sr 
dt p dz po 



dy 




(4) 
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( 5 ) 



dt dz 

Equations (1) through (3) represent u, v, and w momentum 
respectively while (4) and (5) represent continuity and the 
equation of state respectively. They assume an 
incompressible, inviscid linear flow where x, y, z, u, v, and 
w refer to the alongshore, across-shore and vertical 
directions and velocity components respectively, p is 
pressure, is the mean density under the Boussinesq 



approximation , 



ill i 



is the reduced gravity (i.e., assuming 



water column is stratified) , f is the Coriolis parameter, and 
t is time. 

A pressure response form similar to that of Battisti and 
Hickey (1984) 

P = P(z) e'^ gi[iy-o)C] (5) 

is assumed where o) equals wave frequency in radians sec'^. 
Solving for v from equation (2) yields: 

WPo 

Combining equations (3) and (5) produces a second order, 
linear differential equation for w; 
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( 8 ) 



+ M^W = dP ^-kx^illy-<^t] 

dt^ Po 

o 1 dp 

Where N is the Brunt -Vaisala frequency equal to g[-— 

Po 



in 



radians sec'^. Solving for w: 



w = 



20 ) 



dp 



(9) 



p^(iV2-0)2) dz 

The V and w velocities are then eliminated by substituting 
(7) and (9) into (4) producing a second order differential 
equation in P: 



d^P ^ p ^ Q 



dz‘ 



( 10 ) 



(0‘ 



The rigid lid and flat bottom boundary conditions state 
that w = 0 at z = 0 and -H (where H is the bottom depth) . The 
solution to equation (10) takes a sinusoidal form and applying 
the above boundary conditions yields: 



1 ^ _ , 2271 ' ^ 

(02 H 



( 11 ) 



where n is the internal mode number. Finally, solving for : 



0)2 = 






^ 12 
H 



( 12 ) 



This relationship and that for a CSW are used to examine the 
coherent low frequency current fluctuations obseirved during 
Segment 2 of the data set (Chapter V) . 
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II. DATA COLLECTION 



A. CURRENT METER DATA 

Current data were collected at seven different sites from 
Pt. Piedras Blancas to the Farallon Islands. Three separate 
time segments were chosen based on common periods when the 
data overlapped. Tables 1 through 3 summarize these segments. 



TABLE 1. CURRENT METER MOORING LOCATIONS AND DATA 

SOURCES: SEGMENT 1 (17 DEC 89 - 06 APR 90) . 



Mooring 


Bottom 

depth 

(m) 


Meter 

depth 

(m) 


latitude 


Longitude 


Current 

Meter 

Type 


Rotation 
Axis 
C T> 


Data 

Source 


P4 


1090 


350 


35* - 28.3'N 


121" - 44.9'W 


RCM 8 


354 


NPS 


P4 


1090 


500 


35" - 28.3'N 


121" - 44.9'W 


RCM 8 


354 


NPS 1 


P2 


800 


350 


36" - 20.0'N 


122" - 10.2^W 


RCM 8 


344 


NPS 


P2 


800 


500 


36" ‘ 20.0'N 


122" - 10.2'W 


RCM 8 


344 


NPS 


P3 


1800 


350 


36" - 20.0^N 


122" - 27.6'W 


RCM 8 


344 


NPS 


P3 


1800 


500 


36" - 20.0'N 


122" - 27.6^W 


RCM 8 


344 


NPS 


MB1 


650 


500 


36" - 44.3'N 


122" - 02.4'W 


RCM 5 


322 


MBARI 



Since current meter depths from 350 m to 500 m provided the 
most complete and continuous data set, the analysis focused on 
these depths. Although other current meters were deployed, 
only the instruments at the common depths of interest are 
shown. Most of the data were collected with the Aanderaa RCM 
8 current meter with exceptions noted in Tables 1-3. 



13 



Table 2. CURRENT METER MOORING LOCATIONS AND DATA 

SOURCES: SEGMENT 2 (30 AUG 90 - 17 MAR 91) . 



■ = 

Mooring 


Bottom 

depth 

(m) 


Meter 

depth 

Cm) 


Lati tude 


Longitude 


Current 

Meter 

Type 


Rotati on 
axis 
<• T) 


Data 

Source 


P2 


807 


350 


36* - 20. VN 


122* - 10.3'W 


RCM 8 


344 


NPS 


P3 


1800 


350 


36* - 20. VN 


122* - 27.6'W 


RCM 8 


344 


NPS 


MB1 


650 


506 


36* - 44.3'N 


122* - 02.4'W 


RCM 5 


322 


MBARI 


MB2 


1979 


416 


36* - 42.5'N 


122“ - 22.5'W 


RCM 5 


344 


MBARI 



Table 3. CURRENT METER MOORING LOCATIONS AND DATA 

SOURCES: SEGMENT 3 (26 MAY 91 - 18 NOV 91) . 



Mooring 


Bottom 

depth 

(m) 


Meter 

depth 

Cm) 


Latitude 


Longitude 


Current 

Meter 

Type 


Rotation 

axis 

T) 


Data 

Source 


P2(INT) 


807 


400* 


36* - 20. VN 


122* - 10.3'W 


RCM 8 


344 


NPS 


MB1 


650 


478 


36* - 44.3'N 


122* - 02.5'W 


RCM 8 


322 


MBARI 


MB2 


1792 


418 


36* - 40. VN 


122* - 22.5'W 


RCM 8 


344 


MBARI 


I ^ 


1400 


400 


zr - 21.8'N 


123* - 16.0'W 


RCM t 


328 


SAIC 


L_ 


1975 


400 


3r - 39.0'N 


123* • 17.4'W 


RCM 8 


328 


NPS 



* Indicates linear interpolation between 350 m and 500 m 



1. RCM 8 Velocity Measurements 

A large vane assembly (48.5 cm by 50 cm) Icept the 
instriament aligned into the current, A magnetic compass 
inside the recording unit recorded the direction of 
orientation via a potentiometer ring fastened to the lid of 
the compass. Current direction was resolved to 0.35° with an 
accuracy of ± 5° for current speeds from 5 to 100 cm/sec and± 7 .5° 
for speeds from 2.5 to 5 and 100 to 200 cm/sec (Aanderaa, 
1987) . Combined magnetic deviation and variation of each unit 
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was determined on a surveyed test bench at the Naval 
Postgraduate School prior to each deployment. 

A shrouded paddlewheel with magnetic coupling through the 
top end plate measured current speeds in a range from 2 to 250 
cm/sec with accuracies of ± 1 cm/sec or ± 2% of actual speed, 
whichever was greater. The instrument performed a vector 
averaging by sampling the number of rotor revolutions and 
direction every 36 seconds for a 30 minute sampling interval. 
The meter maintained and intermediately stored a running sum 
of North-South and East-West speed components until the data 
sampling interval elapsed. At this time, the resultant 
average speed and its angle were calculated and recorded 
(Aanderaa, 1987) . 

2 . Temperature Measurements 

The RCM 8 used a thermistor with a variable resistance 
(5700 ohm at 0° C and 2001 ohm at 25° C) for temperature 
measurements. The thermistor is connected in series with low 
range (-2.46° C to 21.48° C) , high range (6.08° C to 31.04° C) 
and wide range (-0.34° C to 32.17° C) resistors. Temperature 
values were then determined by various polynomials for 
corresponding temperature ranges based on the measured 
resistances. Temperature readings were resolved to 0.1% of 
the range and accuracy was ±0.06° C. 
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3. RCM 5/RCM 8 Differences 

Current meter moorings MBl, MB2 (Pillsbury, et al . , 
1992) and D deployed the Aanderaa RCM 5 current meter. An 
older version of the RCM 8, the RCM 5 differed in several 
important aspects. First, instead of the solid state data 
storage unit (DSU) used in the RCM 8, the RCM 5 relied on a 
1/4" magnetic tape and tape reader to store and retrieve data. 
Also, the RCM 5 used a Savonius type rotor instead of the 
shrouded paddlewheel used on the RCM 8. The shrouded 
paddlewheel was designed to eliminate the unwanted effect of 
pumping the Savonius rotor through vertical heave in the 
mooring. Additionally, the RCM 8 vane area was redesigned 
from that of the RCM 5 to decrease the response time for 
current direction shifts. Finally, the RCM 5 does not use 
vector averaging. It recorded speed continuously but used a 
snapshot direction at the end of the sampling interval, 
introducing possible errors in current direction measurements. 
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III. ANALYSIS METHODS 



A. NFS DATA 

1. Data Retrieval, Calibration and Editing 

The Aanderaa RCM 8 used a solid state data storage 
unit (DSU) to record current speed and direction, temperature, 
pressure and, when so equipped, conductivity. The DSU had a 
maximum storage capacity of 65,530 ten bit words, which 
corresponds to seven months of data when using a sampling 
interval of 30 minutes. All instruments used a 30 minute 
sampling interval except those at the Monterey Bay moorings 
and mooring D, which sampled data every hour. The RCM 8 
determined sampling time with an internal quartz clock. 
Within all operating temperatures, the clock accuracy was 
±2 sec day'^. 

Upon instrument retrieval, the DSUs were removed and 
downloaded to an IBM compatible personal computer (PC) via the 
DSU reader 2995. The reader converted bit codes to ASCII 
coded format (Aanderaa, 1987) and stored the converted data on 
the PC. Next, a calibration program applied conductivity, 
pressure, speed, time and direction calibrations to convert 
the raw data to engineering units of current speed (cm/s), 
direction (°T) , water temperature (°C) , salinity (PPT) and 
pressure (decibars) . At this point, velocity records were 
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truncated to eliminate unusable data entries i.e., data from 
instrument failures and zero data records obtained before and 
after instrument deployment. The records were plotted and 
visually inspected for obvious outliers. A despiking program 
then checked the data to ensure all parameters were within 
certain tolerance ratios. 

2 . Filtering 

A filtering program (Denbo et al . , 1984) first 
employed a low pass (LP) Cosine- Lanczos filter (half power 
period at 2.9 hours) to remove high frequency noise from the 
data. Signals with periods greater than 3.6 hours passed 
unattenuated and signals with periods shorter than 2 hours 
were removed entirely. The program applied filtering weights 
to successive 25 point segments. A Lagrangian interpolation 
subroutine shifted the start time of the filtered series to an 
even hour. This interpolation scheme employed a third order 
polynomial whose coefficients were based on the desired 
fractional distance between two successive data points. Only 
data with 30 minute sampling intervals were LP filtered. This 
included Farallon (mooring E) , Point Sur and Point Piedras 
Blancas data. 

A low low pass (LLP) filter (half power period at 46.6 
hours) used the output from the LP interpolated filter to form 
six hourly records. This filter, also a Cosine- Lanczos 
type, applied filtering weights to 121 point segments to 
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remove diurnal and shorter period energy as well as any 
inertial period energy. Monterey Bay data had an hourly 
sampling interval already and thus only used LLP filtering. 
The data sets now contained despiked 6 hourly records with 
tidal and inertial influences removed (Figure 3) . 

3 . Bridging 

Gaps were present in various time series due to 
instrument failure or during mooring turnarounds. Gaps of 
four days or less were bridged using a special program 
described below. All discontinuities greater than a few days 
were not filled. The bridge program first determined the 
width of the gap between the two input time series. With 
Lagrangian interpolation, the program then shifted the start 
time of the second series so that it started exactly 6 hours 
after the last time of the gap. A prediction error filter 
(Anderson, 1974) based on data from 25 days on either side of 
the gap determined the data points within the gap. The 
program accomplished this by extending the first time series 
3/4 of the way across the gap forward in time and the second 
time series 3/4 of the way across the gap backwards in time. 
It melded the overlapping data using the prediction filter, 
stored the data separately, and finally combined all data to 
fill the gap. Tests showed that this procedure had no 
detrimental effects on the autospectra. The requirement for 
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Sample Current Meter Filtering: Across -shore 
time series from mooring E 400 m with LP (top) 
and LLP filtering (bottom) . 



20 



APR MAY JUN JIJL 



a continuous time series for spectral and EOF analysis made 
this procedure particularly useful. 

4 . Coordinate Axes Rotation 

The interpretation of geophysical time series is often 
more straightforward in an alongshore vs across -shore rather 
than an east -north coordinate system. Thus, coordinate axes 
rotation facilitated across-shore vs. alongshore analysis. 
This procedure first required primary axis determination 
(Tables 1-3) . A computer program determined the lengths of 
the semi -major axis, semi -minor axis and the angle from the 
positive X axis to the semi-major axis (Kundu, 1976) . The 
program used demeaned u and v to determine the angle from the 
positive X axis, 0: 



0 = 1/2 tan"^ 



N 

— T' (covariance) 

^ bi 



(13) 



where covariance = u^ * v^, u and v are the component speeds 
in the x and y directions respectively and N is the number of 
data points in the record. Based on these results, data were 
rotated locally with the exception of mooring P3 . Because of 
the limited amount of data at P3 and to provide a consistent 
axis in the vicinity of Point Sur, P3 was rotated according to 
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the results from P2 and MB2 . The rotation computed new 
component speeds based on 

= UoCos(4>) + VoSin(4)) 

^rot = 'i^oCos(4)) - u^sin(4>) 

where and are the rotated component speed in the x 

and y direction respectively, u^ and are the initial 

component speeds and <(> is the degrees counterclockwise from 
true north (<() = 90* - |6| from equation (13) above) . The 
rotation process concluded the preliminary data processing 
procedures . 

B . OTHER DATA 

Portions of the data were collected by organizations other 
than the Naval Postgraduate School (NPS) . Specifically, data 
from the Monterey Bay moorings were generously supplied by Dr. 

C. H. Pilskaln of the Monterey Bay Aquarium Research Institute 

(MBARI) . The data were received as calibrated, hourly 

records and were subsequently filtered to form 6-hourly 
records (Pillsbury et al . , 1992) . The data set provided a 

critical piece linking the Farallon Islands and Point Sur 
current meter moorings. Additionally, the Science 

Applications International Corporation (SAIC) worked in 
conjunction with NPS on a project that described the general 
circulation in the vicinity of the Farallon Islands. SAIC 
conducted calibration, maintenance, deployment and data 
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retrieval of all current meters at mooring D. The data at 
mooring D were provided by M. Noble (USGS) , who carried out 
additional quality control, filtering and axis rotation. For 
continuity, all Farallon data (moorings D and E) were rotated 
to 328° T based on these results. 

C. TIME DOMAIN ANALYSIS 

1. Time Series Plots 

Vector velocities, across -shore and alongshore 
velocity components and temperature time series were plotted 
for the common time segments. To facilitate alongshore 
comparisons, stack plots were formed presenting a picture of 
the moorings in a south to north direction; i.e., data from 
the southern most mooring from each segment is on top, the 
northern mooring on the bottom and the others in between. 
These plots proved critical in identifying the various events 
on which to focus additional analysis techniques. 

2. Empirical Orthogonal Fxmction Analysis 

Time domain empirical orthogonal function (EOF) 
analysis (Kundu et al . , 1975; Beardsley et al., 1985; Winant 
et al . , 1987) provided a compact description of the alongshore 
variability of the current structure. The technique found the 
orthogonal eigenvectors and eigenvalues of a covariance matrix 
containing the desired variables. The number of time series 
equalled the number of eigenvalues, with each eigenvalue 
having a corresponding eigenvector or eigenfunction. The sum 
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of eigenvalues equalled the total variance among all 
components of the covariance matrix and thus each eigenvalue 
represented a percentage of the total variance for its 
corresponding mode. Modes were "ranked" according to the 
percent variance contained, i.e., the eigenfunction with the 
greatest variance was called mode 1, the next most mode 2, and 
so on. 

The technique allows one to decompose some parameter 
into a few primary modes representing the essential 
variability. Although meteorologists first employed this 
technique, its applicability to ocean phenomena became quickly 
apparent. Kundu et al . (1975) decomposed current velocity 
fields off the Oregon coast into both dynamical and empirical 
modes . They found that the empirical modes did not depend on 
the accuracy of simplifying dynamical assumptions but only on 
the statistics of the data. They also found that empirical 
modes were discretely orthogonal while dynamical modes were 
not. Beardsley et al., (1985) used both time and frequency 
domains to describe current and temperature variability in the 
Nantucket Shoals Flux Experiment (NFSE79) . Finally, Winant et 
al . , (1987) described moored wind, temperature and current 
observations made during the Coastal Ocean Dynamics 
Experiments (CODE) 1 & 2. In all cases, these researchers 
described the spatial structure of multiple variables in a 
powerful, compact manner. 



24 



3 . EOF Theory 

Following Kundu et al . , (1975) a brief outline of the 

EOF expansion technique is given. For this example, let 
Vj^Cx^) denote demeaned, alongshore velocity at time tj^ (k = 
1, . . . ,K) and at mooring location x^ (i = 1, . . . ,N) . The 

matrix of time series at different mooring locations can be 
decomposed into a modal amplitude and a corresponding 
eigenfunction so that 

N 

= E ^ = 1' ^ 

n=l 

where <1)^ defines the empirical orthogonal eigenfunction and 

Ej^ the modal amplitude for each record in the time series. 
Determining the least squares error between the actual 
velocity and the summation of decomposed velocities shows that 

N 

=X„(|)„(x^), n = l,,,,,N 

i=l 



where 



RiX^.Xj) 



1 

K 



K 



E 

k=l 



v^{x^) v^(x^) 
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represents the real and symmetric covariance matrix, 
corresponds to each eigenvalue and <|)„ obeys the orthogonality 
condition : 

N 

1^1 

where 6^ = 1 for m = n and 0 for m n. The eigenvalues, 

represent the time average energy in each mode i.e., the 
variance and the sum of the eigenvalues both equal the total 
energy. 

A simple MATIiAB function calculated the covariance 
matrix for the desired variables by first removing the time 
mean from each series. The function then found the 
eigenvalues eigenvectors <|)„ and modal amplitudes, as 

described above. 

D. FREQUENCY DOMAIN ANALYSIS 

Spectral analysis employed the Fourier transform to 
provide energy estimates as a function of frequency. A 
modified version of the Matrix Laboratory (MATLAB) function, 
SPECTRUM computed spectral characteristics of across -shore and 
alongshore series. To allow investigation of the maximum 
frequency range without statistical degradation, the data 
segments were divided into five, equal length, overlapping 
pieces. A Hanning window, applied to each piece separately, 
damped the effects of the Gibbs phenomenon resulting from the 
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discontinuity created by the truncation of an infinite series. 
The coefficients of the Hanning window were computed by 

w{n) =0.5(1- cos (2ix-^) ) , n = 1, . . . ,N 

N+1 

where w represented Hanning window weights and N the number of 
data points in each piece. 

Since the piece lengths were not powers of two, the 
modified program applied a mixed-base fast Fourier Transform 
(FFT) to each piece. The transform used by the FFT is: 

N-1 

X{k^l) = ^ x(n+l) 

n=0 

where X defined the transformed time series, x the Hanning 
windowed time series and N the number of records in each piece 
(Little and Shure, 1992) . The program then calculated the 
power spectral density, cross spectral density, complex 
transfer function, coherence and phase functions and 
accumulated results from each piece. The pieces were averaged 
together to provide the final spectral estimates with 
confidence limits. Data were also plotted in variance 
conserving form, obtained by normalizing each spectral 
estimate by the frequency bandwidth and multiplying it by 2.3. 
The results from this and all analysis techniques are 
presented in Chapter IV. 
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IV. RESULTS 



Each time segment presented its own unique results. 
Specifically, currents appeared to vary seasonally between 
each segment. Fluctuations could have been influenced by 
surface seasonal transitions or perhaps deeper current 
variations. This seasonality coupled with distinct, mesoscale 
deviations from the mean flow made it sensible to describe 
each segment separately. 

A. SEGMENT 1 

The first time segment included data from 17 December 19 89 
through 04 April 1990 (110 days) (Figure 2) . Data were 
collected from current meters at P4 off Pt. Piedras Blancas, 
P2 and P3 off Pt. Sur and MBl in Monterey Bay. Instruments at 
350 m and 500 m provided data at all moorings except MBl (500 
m only) . P4 was located approximately 100 km to the south of 
Pt. Sur and 47 km offshore along the 1100 m isobath. Located 
in 650 m of water, MBl was 35 Ion to the north of Pt. Sur, 14 
km offshore and very close to the axis of the Monterey Canyon. 
P2 was positioned 24 km offshore and 13 km from the shelf 
break (200 m isobath) . Figure 4 schematically illustrates 
mooring locations and separation distances. 
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Figure 4. Segment 1 Mooring Separation Distances: 

Alongshore and offshore distances are shown with 
the approximate relative position of the Monterey 
Canyon. The canyon axis is indicated by the 
darker shading. 
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1. Time Domain Analysis 

Basic time domain current statistics for Segment 1 
were calculated for the time period (Tables 4a and 4b) and are 
described for each mooring below. At P4, basic statistics 
showed a mean alongshore (ALS) component speed of 
approximately 1.5 cm sec'^ at both depths with a maximum of 
14.4 cm sec'^ and a minimum of -10.3 cm sec'^ (Table 4a). 
The across -shore (ACS) component mean (-1 cm sec"^) was 
offshore at both depths but is likely within the error of the 
axis rotation. Mean currents were primarily in the alongshore 
direction. The maximum offshore current was -11.9 cm sec'^ 
while the maximum onshore current was 5.9 cm sec'^ (Table 
4b) . From 17 December until 06 January, a vertical shear 
existed at P4 between 350 m and 500 m with poleward flow at 
350 m and equatorward flow at 500 m (Figures 5 and 6) . Flow 
became equatorward at both depths from 06 to 16 January when 
it turned poleward until late winter. Around 08 March, the 
current shifted to equatorward and remained so until the end 
of the segment. The currents looked coherent vertically from 
06 January to the end of the segment. 
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Figure 5. Segment 1 Vector Velocities (350 m) ! Velocity 
vectors are shown in cm sec'^ with poleward 
alongshore, as defined in text, towards the top 
of the plot. Data are from P4 (top) , P2 (middle) 
and P3 (bottom) . 
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Figure 6. Segment 1 Vector Velocities (500 m) : Same as 

Figure 5, except for the 500 m depth level. The 
data from MBl (bottom) are also included in this 
f igure . 
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Table 4a. TIME DOMAIN STATISTICS FOR SEGMENT 1 
(ALONGSHORE COMPONENT) 



! Mooring 


Meter 

depth 

(n) 


Mean 

(cn/s) 


Maxfaui 

(cn/s) 


Miniu 

Ccai/s) 


Standard 

Deviation 

(oa/s) 


Mode 1 

Eigenfmction 
(Eigenvalue = 
177.0 (cm/s)2 
for 66.9% of 
total variance) 


Mode 2 

Eigenfmtion 
(Eigenvalue = 
47.8 (cm/s)2 
for 18.6% of 
total variance) 


P4 


350 


1.5 


14.4 


-10.3 


4.9 


0.22 


0.22 


P4 


500 


1.3 


10.7 


-7.8 


4.0 


0.19 


0.03 


P2 


350 


6.8 


26.1 


-5.9 


7.3 


0.42 


-0.66 


P2 


500 


3.9 


26.0 


-6.7 


5.6 


0.29 


-0.53 


P3 


350 


4.3 


16.8 


-24.4 


9.0 


0.64 


0.41 


P3 


500 


3.4 


11.5 


-17.8 


6.8 


0.49 


0.23 


MB1 


500 


1.2 


8.8 


-7.2 


2.5 


0.04 


0.11 



Table 4b. TIME DOMAIN STATISTICS FOR SEGMENT 1 
(ACROSS -SHORE COMPONENT) 



Mooring 


Meter 

depth 

(n) 


Mean 

(a«/s) 


Maxiaia 

icm/sy 


Niniaui 

(c«/s) 


Standard 

Deviation 

(cm/s) 


Node 1 

Eigenfmtion 
(Eigenvalue = 
45.22 (cm/s) 2 
for 52.8% of 
total variance) 


Mode 2 

Eigenftnction 
( E i genva 1 ue = 
17,51 (cm/sf 
for 20.4% of 
total variance) 


P4 


350 


-1.0 


5.9 


-11.5 


3.8 


0.17 


-0.64 


P4 


500 


-0.7 


3.3 


-10.7 


2.5 


0.17 


-0.37 


P2 


350 


-0.7 


7.4 


-1i.7 


3.4 


-0.12 


0.50 


P2 


500 


0.5 


7.9 


-8.1 


2.7 


-0.12 


0.33 


P3 


350 


-3.4 


9.1 


-16.5 


5.4 


-0.78 


-0.25 


P3 


500 


-2.6 


5.2 


-11.6 


3.8 


-0.54 


-0.16 


MB1 


500 


1.1 


5.2 


-2.0 


1.4 


0.03 


-0.03 



At P2, the mean alongshore current at 350 m for the 
entire segment was 6 . 8 cm sec'^ with a maximum of 28.4 cm sec'^ 
and minimum of -5.9 cm sec'^. The 500 m instrument exhibited 
a similar mean (3,9 cm sec'^) with a maximum of 26 cm sec'^ 



33 



cm sec‘^ and minimum of -6.7 cm sec‘^ (Table 4a) . The across- 
shore component mean was offshore (-0.7 cm sec'^) at 350 m 
with a maximum offshore velocity of -13.7 cm sec"^ and a 
maximum onshore velocity of 7.4 cm sec"^. At 500 m, the mean 
was onshore (0.5 cm sec‘^) (Table 4b) . The small across-shore 
means are again due to the coordinate rotation, which 
statistically minimized the ACS component. 

P2 exhibited much stronger poleward flow than P4 from 
the beginning of the segment until 01 February. On this date, 
the poleward current abruptly diminished, then shifted to 
equatorward for 1 day at 350 m and 10 days at 500 m (Figures 
5 and 6) . As suddenly as the current changed, it shifted back 
to poleward and reached the segment's maximum of 28 cm sec'^ 
at 350 m on 16 February. This short relaxation was not 
observed at the other moorings. The current remained poleward 
until 05 March when the current at both instrument depths 
oscillated with a period of about 10 days for the final month 
of the time segment. The records from 350 and 500 m looked 
vertically coherent the entire time. 

The most consistent poleward flow occurred at P3 . 
With a mean alongshore current of about 4.3 cm sec’^ at 350 
m and 3.4 cm sec"^ at 500 m, the flow remained steadily 
poleward and offshore (across-shore component mean -3.4 and 
-2.6 for 350 m and 500 m respectively) from the beginning of 
the segment until 17 March (Tables 4a and 4b) . On 10 
February, the 350 m current reached a maximum of 17 cm sec'^. 
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This occurred about one week prior to the maximum at P2 . From 
10 February until 17 March, the flow at P3 rapidly shifted to 
equatorward at both depths. Additional data (not shown) show 
that this event exceeded 1000 m depth. During this event, 
current vectors rotated counterclockwise so that the flow was 
directed offshore until 01 April at which time the current 
shifted to onshore until the end of the segment (Figures 5 
and 6) . The maximum equatorward flow exceeded 24 cm sec'^ at 
350 m and 17 cm sec'^ at 500 m. 

With a mean alongshore current of 1.2 cm sec’^ , the 
flow at MBl 500 m demonstrated significantly different 
characteristics throughout the time segment (Table 4a) . The 
poleward currents reached a maximum of 9 cm sec'^ on 03 March 
and a maximum equatorward flow of -7 cm sec'^ on 20 March. 
Weak ( - 5 cm sec'^ ) equatorward bursts appeared with a period 
of about 12 days. Mooring MBl did not display the dominant 
poleward current observed at the other moorings (Figure 6) . 
MBl's location near the axis of the Monterey Canyon seemed to 
isolate the mooring from the general flow pattern offshore 
during this time segment. 

The alongshore and across -shore velocity components 
were separately decomposed into empirical modes using the 
method outlined in Chapter III. Beginning with the alongshore 
component, the first two modes accounted for 85.5% of the 
total variance. The first mode alongshore eigenfunction 
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Figure 7, Eigenfunction Amplitudes for Velocity Modes 1 
and 2 During Time Segment Is ALS (left) and ACS 
(right) eigenfunctions are shown with vectors 
proportionally scaled for visual effect. 
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showed all moorings with the same sign (Figure 7, top left) . 
This represented the strong poleward current seen at moorings 
P2 , P3 and P4 . The second mode alongshore eigenfunction 
indicated an opposing structure between P2 and P3 with P4 
again resembling the flow at P3 (Figure 7, bottom left) . 
This perhaps represents two dynamical regimes (discussed 
later) ; an inshore or upper slope regime (including mooring 
P2) and an offshore or lower slope regime (including moorings 
P3 and P4) . Mooring MBl made only a small contribution to 
these modes, due to its small amplitude and lack of 
correlation with the other moorings. 

For the across -shore component, the first two modes 
accounted for 72.4% of the total variance (Table 4b). The 
first mode eigenfunction indicated opposing flow between Pt. 
Sur moorings (P2 and P3) and the Pt . Piedras Blancas mooring 
(Figure 7, top right) . The second mode eigenfunction however, 
revealed an opposing structure between P2 and P3 with P4 
resembling the flow at P3 (Figure 7, bottom right). This 
again supported the idea that P3 and P4, both located over the 
outer slope were contained within a somehow different 
dynamical regime than P2 and MBl. 

Alongshore modal amplitude and component velocity time 
series were correlated to determine which mooring dominated 
the different eigenmodes (Table 5) . The correlation 
coefficients showed a very high correlation between P3 and 
mode 1 (.94/. 95 at 350 m / 500 m) . The P3 component velocity 
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time series and mode 1 amplitude were plotted and showed 
strong similarities (Figure 8) . This indicated that the 
Table 5. 



P4 350 m 
P4 500 m 
P2 350 m 
P2 500 m 
P3 350 m 
P3 500 m 
MB1 500 m 



SEGMENT 1 CORRELATION COEFFICIENTS FOR FIRST 
AND SECOND MODE EIGENFUNCTIONS WITH TIME SERIES 



Across-shore 


Along 


ishore 


Mode 1 


Mode 2 


Mode 1 


Mode 2 


O.ft 


0.71 


0.59 


0.32 


0.46 


0.62 


0.61 


0.05 


0.24 


0.61 


0.76 


0.62 


0.28 


0.50 


0.69 


0.66 


0.98 


0.19 


0.76 


0.32 


0.96 


0.18 


0.95 


0.24 


0.14 


0.09 


0.23 


0.30 



current at P3 contributed the most to the mode 1 variance. 
Further, mode 2 showed the highest correlation at P2 (not 
shown) suggesting that P2 contributed most to the mode 2 
variance . 

The across -shore modal amplitude time series were also 
computed and correlated with each of the across -shore 
component velocity time series (Table 5) . Consistent with 
the relative magnitudes presented in Figure 7, the first mode 
amplitudes were highly correlated with the P3 time series at 
both 350 m and 500 m (.98 and .96 respectively). The P3 350 
m across -shore component time series was plotted with the mode 
1 amplitude time series clearly showing the current's 
influence on the mode 1 amplitudes (Figure 8) . For mode 2, 
P4 showed the highest correlation (.71 and .62 for 350 m and 
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Figure 8. A Comparison of the Mode 1 Velocity Amplitudes 
with the Best Correlated Individual Time Series 
for Segment 1: The top panel shows the 

alongshore components and the bottom panel shows 
the across -shore components. 

500 m, respectively) wit)i flow at P2 opposing that at P4 and 
slightly less correlated. As expected, the flow at MBl 
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contributed very little to the overall structure in the first 
two modes with correlation coefficients close to zero 
(Table 5) . 

To summarize, time domain analysis described a current 
structure characterized by primarily poleward alongshore flow 
with a small offshore component at P2 and P3 and a small 
onshore component at P4 . A secondary influence caused 
opposing flow between P2 and P3 with P4 generally following 
the P3 pattern. An interesting result showed that although P2 
had a larger mean signature, P3 currents dominated the first 
empirical mode. The significance lies in the EOF analysis 
method where contributions from all moorings were considered 
to produce an overall spatial structure and will be discussed 
in greater detail later. Finally, results from MBl suggested 
that this mooring was secluded from the offshore flow during 
this time segment. 

2 . Frequency Domain Analysis 

Variance conserving spectra were plotted to determine 
the frequencies of maximum energy (Figures 9a through 9d) . 
Since the currents were initially rotated in a manner that 
maximized alongshore variance, the alongshore component 
spectra generally showed the highest energy levels and are 
discussed first. Reported energy levels reflect the 
normalized values to be consistent with the plots. 



40 



VoiiJiiu; Cun^^vlli^ :>pcvtruiu lur F4 35U in 









r. ( c 



Figure 9a. 



Segment 1 Variance Conserving Auto-Spectra for 

350 m: Alongshore and across -shore component 

spectra are shown for P4 (top) and P2 (bottom) . 
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Alongshore Componems Across-shore Components 




Figure 9b. Segment 1 Variance Conserving Auto-Spectra for 
350 m! Alongshore and across -shore components 
are shown for mooring P3 350 m. 
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Figure 9c. Segment 1 Variance Conserving Auto-Spectra for 

500 m: Alongshore and across- shore component 

spectra are shown for P4 (top) and P2 (bottom) . 
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Alongshore Components Across-shore Components 









Figure 9d. Segment 1 Variance Conserving Auto-Spectra 

(500 m) : Alongshore and across -shore component 

spectra are shown for P3 (top) and MBl (bottom) . 
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Alongshore Coinponems Across-shore Components 



At P4 , the 350 m current spectrum showed an energy 
peak of 13.0 cm^ sec'^ at 18 days (Figure 9a, top left) . At 



500 m. 


the 


energy 


decreased 


and 


shifted towards 


shorter 


periods , 


peaking 


at 


8 . 6 cm^ sec"^ 


with a period 


of 


12 days 


(Figure 


9c, 


top 


left) . For 


P2, 


the maximum 


alongshore 



component energy level at 350 m was 28.3 cm^ sec'^ and 
occurred at a period of 18 days . A secondary signal of 15.8 
cm^ sec’^ appeared at a period of 6 days (Figure 9a, bottom 
left) . Similar peaks occurred at 500 m except the energy 
level at 6 days (24.4 cm^ sec’^) actually exceeded the 18 day 
signal of 23.6 cm^ sec*^ (Figure 9c, bottom left). At P3 350 
m, the energy level decreased from that at P2 to a peak level 
of 9.4 cm^ sec'^ at a period of 6 days (Figure 9b, left) . The 
500 m spectrum here also exhibited a maximum energy level at 
a 6 day period but decreased to 5.0 cm^ sec*^ (Figure 9d, top 
left) . The MBl alongshore spectrum (Figure 9d, bottom left) 
displayed a maximum energy level of 6.4 cm^ sec'^ at a period 
of 12 days . 

The across -shore component spectra exhibited energy 
peaks at periods similar to the alongshore components. The P4 
350 m spectrxim peaked at 18 days with an energy level of 14.1 
cm^ sec'^ (Figure 9a, top right) . The energy maximum 
decreased to 4.4 cm^ sec'^ at 500 m and occurred at 12 days 
(Figure 9c, top right) . At mooring P2 , the 350 m spectrum 
showed two peaks again at 12 days and 6 days. The first had 
an energy level of 8.7 cm^ sec*^ while the second was slightly 
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less energetic at 8.0 cm^ sec'^ {Figure 9a, bottom right) . At 
500 m, the energy also peaked at 12 days but decreased to 6.1 
cm^ sec'^. A smaller 3.6 day peak occurred with an energy 
level of 5.0 cm^ sec'^. At P3 , the 350 m spectrum peaked at 
9 days with an energy level of 9.2 cm^ sec'^ (Figure 9b, 
right) . At 500 m, three distinct peaks at 9, 5, and 3.3 days 
occurred with energy levels close to 4.0 cm^ sec'^ at all 
periods (Figure 9d, top right) . The spectrum at MBl 
displayed peaks at 4.5 and 3 days with energy levels around 
2.9 cm^ sec'^ (Figure 9d, bottom right). In general, the 
across -shore component spectra at moorings P4 and P3 displayed 
energy levels comparable to those for the alongshore 
components while P2 and MBl exhibited higher energy level in 
the alongshore component again suggesting two distinct 
dynamical regimes. 

In short, the alongshore component energy levels 
decreased with depth at moorings P4, P2 and P3 with the 
greatest energy at P2 . The energy level at MBl was comparable 
to the 500 m spectra at both P4 and P3 . 

Cross -spectra between moorings showed relatively few 
coherent signals. For the alongshore components (Figures 10a 
and 10b) P2 and P3 were coherent at periods of around 5-6 
days. The signal at P3 led that at P2 by about 40° at 350 m 
and 16° at 500 m. The analysis also showed a coherent 4.5 - 
6 day signal between P4 and P2 at 500 m and a 6 day signal 
between P4 and P2 also at 500 m. P4 led P2 by 26°, 35° and 
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46° for signals of 6 days, 5.1 days and 4.5 days 
respectively. P4 led P3 by 13° for the 6 day signal there. 
The significance of this 6 day fluctuation is not clear since 
it represented such a small fraction of the energy at P4 . 
This may be a spurious coherence resulting from the strong 6 
day peak at mooring P2 . No significant coherence was observed 
at the dominant periods. 

The only across -shore coherent signals were observed 
between P2 and P3 at both 350 m and 500 m. The 350 m data 
(not shown) exhibited coherence at 12 days, 5 days, and 3 
days. At the 12 day and 3 day signals, P2 led P3 by 5° and 
18° respectively while P3 led P2 by 29° at the 5 day signal. 
Coherence values just barely reached the 95% significance 
level at this depth. Although the P2 energy reached a maximum 
at 12 days, the P3 energy at 12 days represented a much 
smaller fraction of the total variance; likewise for the 5 day 
and 3 day signals. Thus the importance of these coherent 
signals is not clear. At 500 m, coherence levels were 
somewhat more significant (Figure 10c) . Specifically, two 
wider coherent bands (9-12 days and 4.5-5 days) were observed 
with an additional peak at 3.6 days. P2 led P3 by about 57° 
for the first band while P3 led P2 by 8° and 48° for the 
second band and by 100° for the 3.6 day signal. These signals 
corresponded with the maximum energy levels observed at P2 50 0 
m (12 and 3.6 days) and those at P3 (9, 5 and 3.3 days). 
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Figure 10a. Segment 1 Coherence and Phase for Alongshore 

Current Components: Coherence and phase are 

shown for P4 with P2 (top) and P4 with P3 
(bottom) . The solid line on the coherence plot 
indicates the 95% significance level. A 
negative phase means that the first series 
leads the second. 
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Figure 10b. Segment 1 Coherence and Phase Between P2 and P3 : 

Coherence and phase are shown for 350 m (top) and 
500 m (bottom) . Other details are the same as in 
Figure lOb. 
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Figure 10c. Segment 1 Coherence and Phase for Across -shore 
Current Components: Details are as in Figure 10a 

except this figure shows the across -shore 
component coherence and phase between P2 500 m 
and P3 500 m only with three peaks from 9-12 
days , 4.5-5 days and 3 . 6 days . 



Therefore, it appeared that the lowest frequency across -shore 
component signal (periods between 9 and 12 days) propagated 
offshore from P2 to P3 while the higher frequencies (periods 
less than 5 days) moved onshore from P3 to P2 for both 
alongshore and across -shore components. 



B. SEGMENT 2 

Segment 2 spanned a time from 30 August 1990 through 17 
March 1991 (200 days) (Figure 2) . Data were collected from 
the Pt . Sur and Monterey Bay moorings only (P2, P3 , MBl, and 
MB2 ) . MB2 was positioned about 30 km seaward of MBl and its 
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position relative to the others is shown schematically in 
Figure 11. The Pt . Sur moorings provided data at 350 m, MB2 
at 416 m and MBl at 500 m. The vertical coherence at mooring 
P2 during other time segments indicated that the currents at 
P2 350 m can reasonably be compared to currents at MB2 416 m. 




Figure 11. Schematic Mooring Locations for Segment 2; 

Distances between moorings are shown as in 
Figure 4. Approximate location of the Monterey 
Canyon is also illustrated. The canyon axis is 
indicated by the darker shading. 



1. Time Domain Analysis 

In general, the mean flow during this time period was 
also directed poleward and offshore (Figure 12, Table 6) . P2 
had a mean alongshore flow of 5.6 cm sec'^ with a maximum of 
23.8 cm sec'^ and a minimum of -10.1 cm sec'^ (Table 6a). 
The across-shore component (mean -0.4 cm sec'^ ) was directed 
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SEGMENT 2; AUG 90 - MAR 91 
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Figure 12. Vector Velocities for Time Segment 2: Velocity 

vectors are shown in cm sec'^ with poleward 
alongshore, as defined in the text, towards the 
top of the plot. Data are from P2 (top), P3 
(second) , MBl (third) , and MB2 (bottom) . 
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Table 6a. TIME DOMAIN STATISTICS FOR SEGMENT 2 
(ALONGSHORE COMPONENT) 



Mooring 


Meter 

depth 

(■) 


Mean 

Ccm/sy 


Maximn 

(OB/S) 


Mininji 

(c«/s) 


Standard 

Deviation 

(cn/s> 


Mode 1 

Eigenfunction 
(Eigenvalue = 
63.22 (cm/s)^ 
for 63.3 % of 
total variance) 


Mode 2 

Eigenfunction 
(Eigenvalue = 
15.83 (cm/s)^ 
for 15.8 % of 
total variance) 


P2 


350 


5.6 


23.8 


-10.1 


6.7 


0.79 


-0.26 


P3 


350 


4.5 


14.5 


-12.0 


4.6 


0.33 


0.93 


MB1 


500 


•0.2 


11.3 


-7.4 


3.1 


-0.13 


-0.10 


MB2 


416 


2.6 


17.7 


-7.8 


5.0 


0.50 


-0.24 



offshore with a maxirtium offshore current of -7.9 cm sec'^ and 
a maximum onshore flow of 5.0 cm sec'^ (Table 6b). For the 
first month, steady poleward flow dominated (Figure 12) . 
Around 28 September, a strongly periodic flow appeared with a 
period between 20 and 30 days. This signal was not clearly 
observed at P3 , located about 25 )an seaward of P2 , but was 
observed at MBl and MB2 . 

The flow at P3 was generally poleward and offshore 
with an alongshore component mean of 4.5 cm sec‘^ with a 
maximum of 14.5 cm sec'^ and a minimum of -12.0 cm sec'^ 
(Table 6a). The across-shore component, with a mean of -2.2 
cm sec‘^ , was again directed offshore. The flow here had a 
strong offshore maximum of -17.5 cm sec'^ around 25 February 
with a maximum onshore current of 7.3 cm sec'^ (Table 6b) . A 
weak (= 7 cm sec'^) equatorward shift appeared around 01 
October and persisted for about 16 days thereafter (Figure 
12) . The current shifted back to poleward around 17 October 



53 



Table 6b. TIME DOMAIN STATISTICS FOR SEGMENT 2 
(ACROSS -SHORE COMPONENT) 



Mooring 


Meter 

depth 

(n) 


Mean 

(OM/S) 


Maxiaui 

icm/s) 


Niniaun 

(ca/s) 


Standard 

Deviation 

(OB/S) 


Mode 1 

Eigenfunction 
(Eigenvalue = 
22.88 (cm/s)2 for 
64.5 % of total 
variance) 


1 

Mode 2 j 

Eigenfmction 
(Eigenvalue = 

5.99 (cm/s)^ for 
16.9% of total 
variance) 


P2 


350 


•0.4 


5.0 


-7.9 


2.4 


-0.27 


-0.12 


P3 


350 


-2.2 


7.3 


-17.5 


4.4 


-0.88 


0.43 


MBl 


500 


1.2 


6.2 


-4.6 


1.5 


0.07 


0.09 


MB2 


416 


0.8 


13.6 


-9.9 


2.9 


0.38 


0.89 



and remained so until 07 March when another equatorward event 
( *= 12 cm sec'^) dominated for the remainder of the segment. 

At MBl, currents were directed equatorward and onshore 
with a mean equatorward flow of 0.2 cm sec'^ and a larger 
onshore mean of 1.2 cm sec"^ (Tables 6a and 6b). Although 
weaker than the other moorings, the stick vector plot (Figure 
12) also showed an oscillatory flow with a period of about 1 
month with phase about opposite (180°) to the currents 
observed at P2 . 

At MB2, alongshore flow was directed poleward with an 
alongshore mean of 2.6 cm sec‘^ , a maximum of 17.7 cm sec'^ 
and a minimum of -7.8 cm sec‘^ (Table 6a) . The across -shore 
flow (0.8 cm sec'^ mean) was slightly onshore exhibiting a 
maximum of 13.6 cm sec"^ and a minimum of -9.9 cm sec"^ 
(Table 6b) . The flow here appeared steadily poleward from 30 
August until 06 November. Beginning around 01 October the 
currents appeared visually coherent with the equatorward 
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bursts described at P2 , but were less energetic. The stick 
vector plots provided initial indications of a coastally 
trapped, poleward propagating disturbance between P2 and the 
Monterey Bay moorings that was less obvious as a modulation of 
the mainly poleward flow at the offshore mooring, P3 . 

EOF results during this time period were again 
consistent with the mean flow statistics. For the alongshore 
components, the first two modes explained 79.1% of the total 
variance. As in Segment 1, the mode 1 eigenfunction showed 
consistency at all moorings except MBl (Figure 13, top left) . 

All first mode alongshore component eigenfunctions had 
the same sign. Together with the stick vector plots, this 
information indicated a tendency toward poleward flow at all 
moorings. The second mode showed a tendency for P2 and the 
Monterey Bay moorings to oppose the flow at P3 (Figure 13, 
bottom left) . This could perhaps be attributed to the low 
frequency fluctuation which was present at P2 and MB2 but not 
at P3 . 

The first two across-shore modes, explaining 80.4% of 
the total variance (Table 6a) , were used to describe the flow 
structure. The mode 1 across-shore eigenfunction described 
offshore flow at the Pt. Sur moorings with P3 supplying the 
largest amount of energy. To the north, an opposing flow 
could be seen at MBl and MB2 (Figure 13, top right) . The 
opposing mode 1 flow between Pt. Sur and MB2 is significant 
and is discussed in greater detail in association with 
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Figure 13. Eigenfunction Amplitudes for Velocity Modes 1 
and 2 During Time Segment 2 : Alongshore 

(left) and across-shore (right) eigenfunctions 
are shown with vectors proportionally scaled 
for visual effect. 



Segment 3. Much like Segment 1, the mode 2 across-shore 
eigenfunction revealed an opposing structure between P2 and 
P3. MB2 followed the structure at P3 (Figure 13, bottom left) 



and appeared as an onshore flow also opposing P2 . The flow at 



56 



MBl contributed more to the total variance than it did during 
segment 1 but on the whole, it played a minor role in the 
overall current structure. 

The mode 1 alongshore amplitude time series showed the 
highest correlation with the P2 alongshore velocity component 
time series (.94) . Thus P2 most nearly resembled the mode 1 
amplitude time series (Figure 14, top), showing P2 as the 
primary source for mode 1 variance. Mode 2 showed the highest 
correlation at P3 (.81) (Table 7). 

Table 7. SEGMENT 2 CORRELATION COEFFICIENTS FOR FIRST AND 



P2 350 m 
P3 350 ID 
MBl 500 m 
M82 416 ID 



SECOND MODE EIGENFUNCTIONS WITH TIME SERIES 



Across-shore 


i 

Alongshore ' 


Mode 1 


Mode 2 


Mode 1 


Mode 2 


0.53 


0.12 


0.96 


0.15 


0.97 


0,24 


0.58 


0.81 


0.21 


0.14 


0.33 


0.13 


0.64 


0.75 


0.79 


0.19 



The across -shore modal amplitude and time series 
correlation showed P3 with the highest correlation (.97) for 
mode 1 while MB2 exhibited the highest correlation for mode 2 
(.75) (Table 7). The P3 across-shore velocity components 
plotted with the mode 1 amplitudes showed this relationship 
very well (Figure 14, bottom) . P2 and MBl did not contribute 
significantly to the across-shore flow spatial structure. 

In general, the primary current structure showed a 
poleward flowing current with a mean slightly smaller than 
that for Segment 1. The second mode alongshore eigenfunction 
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Figure 14. A Comparison of the Mode 1 Velocity Amplitudes 
with the Best Correlated Individual Time Series: 

The top panel shows the P2 alongshore components 
and the bottom panel shows the across -shore 
components . 



revealed consistent flow structures at P2 , MBl and MB2 that 
opposed the P3 structure. The equatorward bursts visible in 
the stick vector plots at these moorings could explain the 
second mode structure. 
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2 . Frequency Domain Analysis 

Variance conserving spectra were calculated for 
alongshore and across -shore velocity components. The 
alongshore component spectra showed energy peaks of 29.2 days 
at all moorings. The energy maximum at P2 (51.2 cm^ sec'^) 
greatly exceeded those at the other moorings (between 9 and 12 
cm^ sec'^) (Figures 15a and 15b) . 

The across -shore component spectra showed significant 
peaks at 19.5 days for P2 and P3 , 29.2 days for MB2 and 3.2 
days for MBl. These peaks were all less than those described 
for the alongshore components (between 4 and 6 cm^ sec'^) 
(Figure I5b) . 

The alongshore components cross -spectra, showed 
coherent signals between several mooring pairs. These 
generally occurred at periods of 58.5, 29.2 and 19.5 days 
(Figures 16a and 16b) . The phase varied with each period and 
between the mooring pairs. Phase diagrams for these three 
periods suggest complex phase relationships (Figure 17) . 
Arrows indicate the direction of phase propagation, e.g. P2 
leads MB2 by 37° at T = 29.2 days. No arrow indicates no 
significant coherence. A marginally significant arrow is 
sometimes shown as a dashed line if the phase was steady and 
consistent with the other pairs. 
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Figure 15a. Segment 2 Variance Conserving Auto-Spectra 

for P2 (top) and P3 (bottom): ALS (left 

panels) and ACS (right panels) at 350 m are 
shown. 
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Figure 15b. 



Segment 2 Variance Conserving Auto -Spectra 
for MBl (top) and MB2 (bottom) : ALS (left 
panels) and ACS (right panels) are shown. 
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Figure 16a. Segment 2 Coherence and Phase for Alongshore 

Current Components: Plots are shown for P2 

with MBl (top) and P2 with MB2 (bottom) . 
Other details are as in Figure 10a. 
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Figure 16b. Segment 2 Coherence and Phase for Alongshore 

Current Components: Plots are shown for MBl 

with MB2 (top) and P3 with MB2 . Other details 
are as in Figure 10a. 
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The 29.2 day signal (Figure 17, middle) represented 
the largest energy peak. At this period, P2 led P3 by 48° and 
MB2 by 38° consistent with a poleward propagating disturbance. 
Another interesting observation could be made by inspecting 
the phase between P2, MB2 and MBl. P2 led MB2 by around 40° 
and MB2 led MBl by about 135°. Also, consistent with the 
stick vector plots, MBl appeared to be 180° out of phase with 
P2 . Adding the phase from P2 to MB2 with that from MB2 to 
MBl, yields a result consistent with the phase between P2 and 
MBl and adds confidence to the results. 

The 19.5 day signal showed some characteristics that 
were consistent with the other two periods. Specifically, P2 
led both MB2 and P3 by 42° and 23°, respectively. In contrast 
to the other signals, MB2 led P3 by 24° and MBl led P2 by 
170°. The energy at 19.5 days was less than that of 29.2 days 
but significantly higher than the 58.5 day signal. 

For the 58.5 day signal, P2 led MB2 by 42°, MBl by 
170°, and MB2 by 23°. P3 was coherent with MBl and led by 
136° and MB2 led MBl by 132°, much like the 29.2 day signal. 
In general, energy seemed to propagate in a poleward direction 
but the significance of this picture is probably small since 
the energy level at 58.5 days was the lowest of all coherent 
periods . 

By knowing the phase and mooring separation distances, 
alongshore wavelengths between P2 and MB2 and across -shore 
wavelengths between P2 and P3 could be calculated. The 
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T = 58 Days 



T = 29.2 Daya 



1 32 





Figure 17. Segment 2 Phase Dlagramt Arrows indicate the 
direction of phase propagation between moorings. 
The phase in degrees is shown for each period of 
significant coherence. 



results for the most significant energy levels (29.2 and 19.5 
days) are summarized in Table 8. The phase error of cross - 
spectral analysis are also presented here. It is important to 
note the phase errors since this could introduce a significant 
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source of error in the wavelength calculations, which will be 



compared to theory in Chapter V. 



Table 8. OBSERVED PHASE AND WAVELENGTHS FOR LOW FREQUENCY 



T = 29.2 days 
T = 29.2 days 
T = 29.2 days 
T = 19.5 days 
T = 19.5 days 



OSCILLATION (SEGMENT 2) : A positive phase speed 
indicates that the first mooring leads the second. 



Mooring 
Pai r 


Phase 

(0) 


Phase 

(days) 


Phase 

error 

(0) 


Phase 
speed 
(km/ day) 


Observed 

wavelength 

(km) 


Separation 

Distance 

(km) 


P2 & MB2 


38 


3.1 


49 


+13 


389.5 

(ALS) 


41 


P2 & P3 


48 


3.9 


34 


+7 


192.1 

(ACS) 


26 


MB1 & M82 


137 


11.1 


46 


-3 


78.8 

(ACS) 


30 


P2 & MB2 


42 


2.3 


26 


+ 18 


351.0 

(ALS) 


41 


P2 & P3 


75 


1.2 


41 


+6 


122.3 

(ACS) 


26 



Observed wavelengths were calculated by multiplying 
the separation distances by 360° and dividing by the phase 
difference. This procedure produced alongshore wavelengths of 
390 km for the 29.2 day signal and 351 km for the 19.5 day 
signal. The across -shore wavelengths were 192 km and 122 km 
for the 29.2 and 19.5 day signals, respectively at Pt. Sur. 
Between the Monterey Bay moorings, an across -shore wavelength 
of 78.8 km was calculated for the 29.2 day signal but with a 
phase speed in the opposite direction from off Pt. Sur (Table 
8) . This procedure calculated the largest possible wavelength 
based on the phase differences. The phase speed was 
calculated by dividing the alongshore separation (41.1 km 
between P2 and MB2) and the across-shore separation (26 km 
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between P2 and P3 , 30 km between MBl and MB2) by the time 

lags . 

Cross spectra calculations showed little coherence 
between the across -shore velocity components. The only 
significant coherence occurred between moorings P2 and MB2 at 
a period of 19.5 days (not shown) . MB2 led P2 by about 175° 
so the two moorings were generally out of phase in the across - 
shore component. A similar situation existed in Segment 3, 
where a mesoscale meandering feature dominated the currents at 
these moorings. 

C . SEGMENT 3 

Segment 3 lasted from 26 May 1991 until 18 November 1991 
(177 days) (Figure 2) . In addition to the Pt . Sur and 
Monterey Bay moorings, this time segment included data from 
moorings D and E near the Farallon Islands. Mooring D was 
located about 110 )on to the north of MB2 and mooring E was 
positioned 32 km. to the north of D (Figure 18) . Currents at 
P2 were obtained at depths of 350 m and 500 m but the other 
moorings only had data at 400 m. For purposes of EOF and 
spectral analyses, the currents at P2 were linearly 
interpolated to a depth of 400 m to approximately match the 
others (Figure 19) . 

1. Time Domain Analysis 

Time domain statistics again showed a poleward 
alongshore component mean at all moorings. At P2 , the mean 
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Figure 18. Mooring Separation Distances for Segment 3: 

Details are as in Figure 4. 

alongshore current was 6.3 cm sec'^ with a maximum of 33.3 cm 
sec'^ and a minimum of -19.0 cm sec"^ (Table 9a). The 
across -shore statistics (Table 9b) showed a mean offshore flow 
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Figure 19. 



Segment 3 Vector Velocitlees Velocity vectors 
are shown as in Figure 5. Data are from P2 
(interpolated to 400 m; top) , MBl (second from 
the top) , MB2 (third from the top) , D (second 
from the bottom) and E (bottom) . 
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of -1.0 cm sec'^ with a maximum onshore flow of 9.8 cm sec'^ 
and a maximum offshore flow of -18.6 cm sec"^. 

The statistics at MB2 were similar showing an 
alongshore mean of 4.2 cm sec'^ with a maximum of 23.2 cm 
sec*^ and a minimum of -10.2 cm sec'^ (Table 9a) . The 
across -shore component had a mean of 0.9 cm sec*^ with a 
maximum of 17.9 cm sec"^ and a minimum of -10.9 cm sec'^ 
(Table 9b) . 



Table 9a. TIME DOMAIN STATISTICS FOR SEGMENT 3 
(ALONGSHORE COMPONENT) 



Mooring 


Meter 

depth 

(■) 


Mean 

(an/s) 


Maxiaui 

Ccm/sy 


Minimi 

<Oi/S> 


Standard 

Deviation 

(ca/s) 


Mode 1 

Eigenfunction 

(Eigenvalue = 
193.80 (cm/s)2 
for 64.7% of 
total variance) 


Mode 2 

Eigenfunction 
(Eigenvalue = 
65.18 (cm/s)^ 
for 21.7% of 1 

total variance) | 


P2 


400* 


6.27 


33.3 


-19.0 


10.73 


0.70 


-0.47 


MB1 


500 


-0.56 


8.4 


-8.7 


2.78 


0.02 


-0.07 


MB2 


416 


4.19 


23.2 


-10.4 


6.19 


0.30 


-0.38 


D 


400 


5.99 


23.6 


-6.3 


6.42 


0.36 


0.07 


E 


400 


3.64 


24.6 


-21.8 


9.84 


0.53 


0.79 



* velocities interpolated between 350 m and 500 m 



From the beginning of the time segment until 27 June 
the flow was poleward at P2 400 m, with the maximum over the 
entire segment (33.3 cm sec'^) occurring around 02 June. On 
27 June, the flow abruptly shifted to equatorward for four 
days and then reversed back to poleward until early August. 

From around 07 August to 07 September a relatively 
strong (10-20 cm sec'^) equatorward event dominated the 
currents at both P2 and MB2 . The across-shore components at 
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these two moorings exhibited a remarkably strong negative 
correlation. Examination of the across -shore components of P2 
and MB2 alone (Figure 20) , clearly showed onshore flow at MB2 
opposing offshore flow at P2 suggesting the presence of a 
mesoscale, rotating feature between the two moorings. 



Table 9b. TIME DOMAIN STATISTICS FOR SEGMENT 3 
(ACROSS -SHORE COMPONENT) 



Mooring 


Meter 

depth 

(n) 


Mean 

(cm/s) 


Maximum 

(cm/s) 


Miniiui 

(cm/s) 


Standard 

Deviation 

(cm/s) 


Mode 1 

Eigenfinction 
(Eigenvalue = 
31.96 (cm/s)2 
for 41.3% of 
total variance) 


Mode 2 

Eigenfunction 
(Eigenvalue = 
19.45 (cm/s)^ 
for 25.1% of 
total variance) 


P2 


400* 


-1.0 


9.8 


-18.6 


4.9 


-0.06 


t 

0.16 


MB1 


500 


1.2 


6.4 


-3.6 


1.3 


-0.06 


0.03 


MB2 


416 


0.9 


17.9 


-10.9 


4.0 


0.56 


-0.01 


D 


400 


0.8 


19.3 


-12.9 


4.2 


-0.17 


-0.67 


E 


400 


-0.6 


10.4 


-19.0 


4.2 


-0.04 


-0.73 



* velocities interpolated between 350 m and 500 m 



P2 currents in early August prior to the event were 
poleward but rotated clockwise slowly until 03 August. The 
current then shifted abruptly to equatorward and onshore 
briefly. Finally, it turned offshore on 07 August to begin 
the event. The current maintained an offshore, equatorward 
direction at P2 until the event ceased around 07 September 
(Figure 21) . Following the event, the current shifted to 
poleward and onshore for approximately one month until an 
equatorward burst of about -10 cm sec'^ interrupted the flow 
on 28 September (Figure 19) . A similar deviation occurred on 
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28 October. These equatorward signals appeared at about the 
same period as those observed in Segment 2 (30 days) . 

Currents at MB2, directed poleward and onshore during 
the beginning of the event, did not appear to feel the 
feature's influence until 13 August. On this day the current 

II (’nm|M)ncn( Vol(K:iiics l(»r l'2 and MU2 
MII2 116 m 




*20 



Figure 20. 
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Across -shore Component Velocities at P2 and 
MB2: A mesoscale event commenced around 07 

August, characterized by negative correlation 
for the next 25 days. 



turned to equatorward and onshore and remained so until 07 
September (Figure 21) . Flow also shifted to poleward after 
the event here and followed a similar pattern as that 
described for P2 ; i.e. a poleward relaxation to a near 0 cm 
sec'^ current on 28 September followed by an equatorward 
shift (-5 cm sec'^) 30 days later. The flow at both P2 and 
MB2 during October and November appeared remarlcably similar to 
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chat during the same months, one year earlier in Segment 2 
( Figure 12 ) . 

In general, poleward flow dominated from May until 
early August at both Farallon moorings. At mooring E though, 
the flow shifted to equatorward (max. -18 cm sec'^ ) on 22 
June and remained so for about 9 days thereafter. This 
occurred about five days prior to the equatorward shift 
described at P2, approximately 183 km to the south. However, 
moorings MB2 and D located between P2 and E, did not show such 
a shift so a connection between these equatorward 
interruptions does not seem likely. The shift at mooring E 
was probably due to a cold filament observed off Pt . Reyes, 
north of San Francisco Bay during this time segment (Ramp et 
al . , 1992). Around 21 August, approximately two weeks after 
the start of the event at P2 and MB2 , flow at moorings D and 
E weakened significantly and shifted from poleward to 
equatorward. This raised the question of whether or not a 
mesoscale event occurring some 140 km to the south could 
manifest itself into the flow pattern at D and E. With the 
limited data available, a connection between these two events 
seemed speculative but it does present an interesting 
possibility. 

At mooring D, the flow shifted back Co poleward around 
27 September and remained weak (< 6 cm sec'^ ) for the 
remainder of the segment (Figure 19) . At mooring E, a strong 
equatorward reversal occurred around 23 September, which was 
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Figure 21. Segment 3 Vector Velocities During Mesoscale 
Event! Velocity vectors are shown 25 July 
through 04 September with details as in 
Figure 5. Data are from MB2 (top) and P2 
(bottom) . 
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not observed at mooring D. For the next month, the current 
at E grew steadily to a maximum of -18 cm sec"^. It then 
shifted back to poleward on 21 October and remained fairly 
weak and much more variable (Figure 19) . 

In the empirical mode decomposition, the first two 
alongshore modes accounted for 66.4% of the variance. The 
mode 1 alongshore eigenfunction (Figure 22, top left) showed 
the flow at all moorings to have the same sign. As discussed 
previously, this could be interpreted as the poleward 
undercurrent seen in the stick vector plots. The mode 2 
alongshore eigenfunction showed a divergent flow between the 
Farallon moorings and those to the south (Figure 22, bottom 
left) . This could perhaps be due to the rotating feature 
observed at P2 and MB2 , which impacted the currents there but 
not farther north at moorings D and E. 

As with the other segments, the mode 1 across -shore 
eigenfunction at P2 and MB2 described an opposing flow 
structure (Figure 22, top right) . Moorings D and E were not 
strongly correlated with the Pt. Sur and Monterey Bay moorings 
indicating that the mesoscale event at these southern 
locations dominated the mode 1 across -shore component 
structure. The mesoscale feature described earlier appeared 
to influence the mode 2 across -shore current structure as 
well. The MB2 mode 2 eigenfunction had a sign opposite to P2 
(Table 9b) but the magnitude (-0.01) was much smaller than 
that at P2 (0.16) and thus does not show in Figure 22. 
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Modal amplitudes were again correlated with alongshore 
component time series (Table 10) . P2 showed the strongest 
correlation with the mode 1 amplitudes having a correlation 
coefficient of 0.91. The P2 alongshore component plot and 
modal amplitude time series showed that this mooring 






Figure 22. Eigenfunction Amplitudes for Velocity Modes 1 
and 2 During Time Segment 3: Alongshore (left) 

and across-shore (right) eigenfunctions are 
shown with vectors proportionally scaled for 
visual effect. The across-shore modal 
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Table 10. 



SEGMENT 3 CORRELATION COEFFICIENTS FOR FIRST 
AND SECOND MODE EIGENFUNCTIONS WITH TIME SERIES 



P2 400 m 
MB1 500 m 
MB2 418 m 
D 400 m 
E 400 m 



Cross-shore 


Alongshore 


Mode 1 


Mode 2 


Mode 1 


Mode 2 


1 0.91 


0.35 


0.92 


1 

0.15 ' 


0.10 


0.20 


0.20 


0.('5 


0.67 


0.50 


0.80 


0.01 


0.78 


0.09 


0.23 


0.70 


0.75 


0.65 


0.05 


0.76 



contributed most of all the moorings to the mode 1 variance 
(Figure 23) . Mooring E was most nearly correlated with the 
second mode alongshore eigenfunction (Table 10) . P2 again 
showed the strongest correlation (0.92 coefficient) and the 
relationship between the two is illustrated in Figure 23. 
Much like the alongshore components, mooring E was best 
correlated with the mode 2 eigenfunction. 

2 . Freqfuency Domain Analysis 

The variance conserving spectra (Figures 24a through 
24 c) for alongshore components showed that Moorings D, E and 
P2 peaked at 29.2 days, similar to Segment 2. MB2 experienced 
a maximum energy level at a shorter period of 14.6 days . 
Mooring E contained the most energy, 102.9 cm^ sec'^. P2 had 
less energy with a 68.4 cm^ sec’^ peak. MBl showed a 19.5 day 
peak with an energy level of 11. 0 cm^ sec'^. D and MB2 peaked 
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Figure 23. A Comparison of the Mode 1 Velocity Amplitudes 
with the Best Correlated Individual Time Series: 

The top panel shows the alongshore components 
and the bottom panel shows the across -shore 
components . 
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Figure 24a. 



Segment 3 Variance Conserving Auto-Spectra for 
P2 (top) and MBl (bottom) . 



79 



Aloiiushore Coinponems Across-shore Components 



Vanancc 0*nsi;rvtnc Sncurum inr MB2 4 18 in Vanance Cuiwcrvmg Sprciruin lor MB2 418 in 












Figure 24b. Segment 3 Variance Conserving Auto-Spectra for 
MB2 (top) and D (bottom) . 
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Alongshore Componenis Across-shore Components 








Figure 24c. Variance Conserving Auto-Spectra for Segment 3 

Mooring E 
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Alongshore Coinponenis Across-shore Components 



respectively 



with energy levels of 9.0 and 17.6 cm^ sec 
(Figures 24a through 24c) . 

Across -shore component energy peaks occurred at 
periods ranging from 5.8 days (moorings D and MBl) , 7.3 days 
(moorings E and MB2) and 29.2 days at P2 (Figures 24a through 
24c) . The peak at P2 contained the maximum energy of all 
moorings (39.0 cm^ sec'^) (Figures 24a through 24c). The 
pervasiveness of the 29.2 day peak is interesting, but missing 
this time at MB2 in both the alongshore and across -shore 
components . 




Figure 25. Segment 3 Coherence and Phase for Across -shore 
Components for P2 and MB2 : The solid line on 

the coherence plot indicates the 95% 
significance level. The two mooring locations 
are coherent but 180° out of phase. 
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Cross spectra showed the only coherent alongshore 
component signal between moorings D and E at periods from 9.7 
to 8.4 days. The energy levels at these periods were 
significantly smaller than at longer periods; thus the 
importance of this coherence is not clear. Significant 
coherence occurred between the P2 and MB2 across -shore 
components at periods from 58 to 19.5 days (Figure 25) . These 
signals were =180° out of phase supporting the opposing flow 
structure resulting from the presence of a rotating mesoscale 
feature. No other coherent signals were observed between any 
other mooring pairs. 
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V. DISCUSSION 



A. CALIFORNIA UNDERCURRENT 

Time series plots and EOF analysis consistently 
demonstrated a predominantly poleward flow during the three 
time segments at all moorings. These results were generally 
consistent with the flow of the California Undercurrent (CUC) 
as described by Chelton (1984), Wickham et . al (1987) and 
Robson (1990) . Chelton (1984) computed seasonal geostrophic 
velocities and found the maximum poleward current occurred 
during December off Point Sur. Segments 1 and 2 time series, 
however, showed current maxima during February. Chelton based 
his calculations on long term averages of all available data 
from the California Cooperative Oceanic Fisheries 
Investigations (CALCOFI) . The reliability varied from month 
to month and from location to location and could possibly 
explain the difference between Chelton' s findings and those of 
Segments 1 and 2. Also, Chelton considered the entire water 
column to include currents extending to the surface while 
these results examined the 350 - 500 m depths only. The 
important point is that definite poleward maxima existed 
during winter in both Segments. 

The CUC signature during Segment 3 agreed well with 
Wickham et. al (1987) who found predominantly annual 



84 



variability at Point Sur with a maximum poleward flow in May 
and June. Specifically, Segment 3 time series clearly showed 
a maximum at all locations to include Point Sur, Monterey Bay, 
and the Farallon Islands moorings during May 1991. During the 
first two segments, the flow at P3 provided the most 
consistent poleward flow and did not exhibit the periodic 
reversals to equatorward flow seen at P2 . 

B . COASTAL TRAPPING 

Although the poleward flowing CUC dominated the mean 
currents, a low frequency oscillation was superimposed on the 
poleward undercurrent during Segment 2 (October 1990 to 
January 1991) . This deviation was observed more clearly at P2 
than at P3 indicating the presence of a coastally trapped 
phenomenon. If a wavelike feature is assumed, the coherence 
and phase information presented earlier can be used to 
calculate phase speeds and wavelengths for the different 
energetic peaks in the spectrum (Table 8) . These can be 
compared with the dispersion relations for different kinds of 
simple coastally trapped waves to see if they provide a 
reasonable explanation for the frequencies, wavelengths and 
phase speeds observed. 

1. Rossby radius of deformation 

Significant evidence supported the importance of 
coastal - trapping . First, calculation of the internal Rossby 
Radius of deformation revealed that P2 lied close to one 
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Rossby radius of the shelf break while P3 was located outside 
of this radius. The internal Rossby radius of deformation 
represents the horizontal scale at which rotation effects 
become as important as the buoyancy effects (Gill 1982) . An 
expression for this scale can be derived following the 
procedure outlined in Gill (1982) : 



where 

= M. 

gn-n. 




= represents the wave speed of the n*^*^ 

is the equivalent depth. An expression 



( 14 ) 

mode and 
for the 



Rossby radius was determined by substituting the equivalent 
depth into the equation for and then substituting into 

(14) to get: 



N H 






mtf 



( 15 ) 



CTD casts conducted near P2 provided Brunt- Vaisala (N^) 
profiles for October, November, December 1990 and February, 
1991. Using the standard equation: 



p OZ 



a depth averaged value of 5.50 x 10'^ rad sec"^ was found. 
Using the current meter mooring depth at P2 (800 m) for H and 
the Coriolis parameter of 8.65 x 10'^ rad sec'^ produced a 
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first mode (n=l) Rossby radius of 16.2 km, and explained why 
coastal trapping influenced the flow at P2 more strongly than 
it did at P3 (Figure 26) . 

DISTANCE OFFSHORE (KM) 

ROSSBY RADIUS 
(16 KM ) 




2000 



Figure 26. Schematic Cross Section Between P2 and P3 : 

Diagram shows the bottom slope and the Rossby 
radius of deformation, drawn to scale at the top 
of the plot. 
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2. Critical Latitude for Coastal Trapping 

Additional support for coastal- trapping in the 
vicinity of Point Sur came from a study by Clarke and Shi 
(1991) . A simple formula produced the critical frequency at 
which waves become trapped at a given latitude: 

C*)<- = COS0 (16) 



where P is the latitudinal change in planetary vorticity 
(2 X 10'^^ m'^ sec'^) , Cj^ = the phase speed for the n*-*^ internal 
mode, f^ = the Coriolis parameter (8.64 x 10'^ sec'^) and 0 is 
the angle between the coastline and true north determined from 
a chart to be 35°. Using the expression: 



with N = 5 .5 X 10'^ rad sec'^ and H = 800 m, gives a value of 
C]^ = 1.40 m sec'^. Substituting into (16) yields a frequency 

of 1.33 X 10"^ sec'^ or a period (T=-^) of 548 days. Waves 

0 ) 

with periods greater than this will leak offshore in the form 
of a planetary Rossby wave while smaller periods are coastally 
trapped. Both observed periods (29.2 and 19.5 days) were well 
within this maximum trapping period offering additional 
evidence to strengthen the argument that the low frequency 
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oscillations observed off Point Sur were trapped near the 
coast . 

3 . Poleward Propagation 

Alongshore component phase observations from 
Segment 2 show that the oscillation propagated poleward at all 
periods (Figures 16 and 17) . This supports the theory- 
presented in Chapter I that coastally trapped phenomena 
propagate in a direction such that the coastal boundary is on 
the right in the northern hemisphere. In all cases except the 
19.5 day signal, the equatorward reversals propagated from P2 
towards the north (Figure 17) . This evidence along with the 
Rossby radius information and critical trapping latitude 
prompted a more detailed analysis of possible coastal trapping 
mechanisms. The analysis focused on alongshore component 
phase observations only, since the majority of the variance 
was contained here. 

C. Possible Coastal Trapping Processes 

1. Internal Kelvin Waves 

As discussed in Chapter I, the presence of an IKW 
first required the assumption of a flat bottom with a vertical 
boundary at the coast. The solution to the momentum equations 
based on these assiamptions produced two waves: one that 

exponentially decayed from the coast seaward and the other 
which grew exponentially from the coast seaward. The former 
represents the only physical solution and propagates poleward 
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along an eastern boundary. The dispersion relationship of 
Equation (12) facilitated a comparison between the observed 
wavelengths and theoretical wavelengths. Only the periods of 
maximum energy (29.2 and 19.5 days) were considered. Since o) 
and H were well known compared with 1, the dispersion relation 
was rearranged to produce a solution for the alongshore wave 
number for each internal mode (n) : 



J2 



(i^o» 

H 



{N 



2 _ 



o>2) 



( 17 ) 



By substituting known quantities into (17) , a theoretical 
value for the alongshore wave numbers and hence wavelengths 
for the two periods of interest could be obtained. In (17) : 

= 29.2 days, o) = = 2.49 x 10‘® rad sec'^ (T 2 = 19.5 



days, ( 1)2 = 3.70 x 10'® rad sec'^) , H = the P2 mooring depth 

(800 m) and N = 5.5 x 10’^ rad sec’^, i.e., the same value 
used in the Rossby radius calculation. The results (Table 11) 
show that IKWs of these periods are expected to have much 
larger wavelengths than those observed. The wavelengths 
decrease with increasing n, but modes higher than 1 or 2 have 
not been observed in the ocean so this seems an unlikely 
explanation. Our N is, if anything an underestimate, and 
larger H makes the theoretical wavelength larger so this 
cannot reduce the discrepancy. 
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Table 11. OBSERVED WAVELENGTHS COMPARED WITH INTERNAL 



T, = 29.2 Days 
T2 = 19.5 Days 

Since N was well observed, we conclude that pure IKWs cannot 
explain the observed phase relationships during Segment 2. 

2. Barotropic Continental Shelf Waves 

These waves provide the opposite "end point" for 
comparison with the calculations. As in the IKW procedure, 
rearranging the CSW dispersion relationship from Chapter I 
provided a means to calculate a theoretical alongshore 
wavenumber and wavelength for a CSW: 

+ ;c2 +12 (18) 

0 ) 

where f is the Coriolis parameter (8.64 x 10^ rad sec' for 
36°-20' N) , k is the observed across-shore wave number, X is 

the bottom slope parameter (-— — and o) is the observed 

2/3 QX 

wave frequency. Using a depth of 800 m (P2 mooring depth) , 
a dh of 1800 m, a dx of 52 km (distance to P3) , X = 

2.16 X 10 ’^ m' . Alongshore wavenumber calculations produced 
a wavelength of 28 km for the 29.2 day signal and 27 km for 
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the 19.5 day signal (Table 11). Assuming that less than one 
complete wavelength fit between mooring locations, these waves 
are too short to explain the observations. The preliminary 
conclusion is that the waves are either 1) hybrid waves where 
both bottom slope and stratification are important or 2) waves 
which are scattered/reflected by the Monterey Canyon, which 
renders conventional dispersion relations useless. 

3. Limitations of CTW Theory 

The differences in Table 11 can be explained by 
examining Che governing assumptions behind each dispersion 
relationship. First, IKW theory assumed a flat bottom with a 
vertical wall. This assumption is not valid since the 
Monterey Canyon passes directly between P2 and MB2 (Figures 1 
and 11) . On the other hand, CSW theory assumed homogeneity 
with a sloping bottom. Hydrographic data in this region shows 
vertical density stratification (Robson, 1990) and weakens Che 
homogeneous assumption. Therefore, it seems Chat the observed 
oscillations are produced by a complex interaction of both 
processes . 

Finally, while the wavelength calculations for IKW is 
non- ambiguous , it is possible in the case of shelf waves Chat 
more Chan one wavelength "fits" between moorings P2 and MB2 . 
If Che actual phase difference, for instance, was 360° + 38° 
= 398°, rather than 38°, then a wavelength of 37 Jon is 
implied, in closer agreement with the observed wavelength. 
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This idea is speculative, but does represent an alternate 
solution to the problem. 

4. Additional Evidence Supporting Coastal Trapping 

Vertical coherence and phase were calculated at P2 
from May 19 89 through October 1990 (17 months) . Current meter 
data at 100 m, 350 m and 500 m showed an upward vertical phase 
propagation consistent with results of Romea and Allen (1983) 
suggesting an IKW presence. Romea and Allen (1983) focused on 
low frequency behavior with periods similar to those observed 
during Segment 2 (T = 20 - 40 days) . They employed a model 
that examined "the response of a rotating, stratified, f-plane 
ocean with a rigid lid, forced by an alongshore wind stress at 
the coast." They calculated the response for three different 
types of wind stresses: winds acting for a brief time at a 

single point, winds acting for a longer time at a single point 
and more generally, winds acting over a larger area for a long 
time. Their model showed that bottom reflections could be 
ignored if (1) the bottom is highly dissipative, (2) if the 
waves have sufficiently low frequencies and travel a great 
distance before hitting the bottom or (3) if the waves are 
internally dissipated prior to reaching the bottom. They also 
showed that a poleward travelling wind forced a poleward 
propagating component in the form of a coastally trapped IKW 
that had a phase velocity directed upward and a group velocity 
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directed downward i.e., consistent with the results observed 
at P2. 

Several other models have addressed CTW forcing 
mechanisms. Chapman (1987) applied wind forced, long coastal- 
trapped wave theory in an attempt to hindcast quantities 
observed during the Coastal Ocean Dynamics Experiment (CODE) 
which took place in the vicinity of Point Arena during 1981 
and 1982. He found that bottom pressure and alongshore 
current were forced by winds between Pigeon Point, (200 1cm 
south of the CODE site) and the CODE site itself. Battisti 
and Hickey (1984) also found that sub surface pressure (SSP) 
response off Washington and Oregon resulted from wind forcing 
between San Francisco and Cape Mendocino over the entire low 
frequency band for periods greater than 5 days. They found 
that local wind forcing was negligible during the summers of 
1972 and 1978, and thus remote forcing was responsible for 
fluctuations in SSP observed then. Finally, Lopez and Clarke 
(1989) considered a solution forced by combined local and 
remote winds where the local solution disregarded the 
alongshore pressure gradient but allowed local winds and the 
remote solution assumed local winds were zero and the 
alongshore pressure gradient Py 0. They considered periods 
between 22 and 45 days. Although these models focused on 
continental shelf phenomena, the results consistently show 
that remote wind forcing, south of the current observation 



94 



resulted in coastal trapping at periods in general agreement 
with those observed during Segment 2 . 

Vertical coherence and phase were also calculated at 
P3 with data from October 1990 through May 1991 (8 months) . 

During this time period current meters were located at depths 
of 100 m, 350 m, 500 m and 1000 m. In contrast to the results 
at P2 , a downward vertical phase propagation was observed for 
all current meter combinations. This indicates that P3 was 
more affected by surface phenomena such as meanders and 
eddies, than by coastal trapping. 

As discussed in the previous chapter, other 
significant dynamical differences between P2 and P3 were 
observed during the first two time segments. To review, the 
current at P3 during Segment 1 showed a persistent poleward 
flow interrupted only by a strong equatorward event toward the 
end of the segment. This event appeared at the surface and 
extended down to at least 1000 m (unpublished data) . 
Additionally, the equatorward reversals observed at P2 and the 
Monterey Bay moorings during Segment 2 were less apparent at 
P3, i.e., currents at P3 appeared to be influenced less by the 
coastal boundary. This makes sense in light of P3's location 
well outside of the Rossby radius. With this evidence, we 
conclude that data from P2 and P3 were representative of two 
dynamically unique environments: an upper slope or near shore 

regime (P2) and a lower slope or offshore regime (P3) . 
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5. Further Possible Analysis 

A more rigorous analysis of this problem would require 
a dedicated numerical model that considers both irregular 
bathymetry and density stratification. Wilkin and Chapman 
(1990) used such a model to study freely-propagating CTWs. 
Their model attempted to simulate scattering of CTWs at a site 
on the East Coast of Australia. They found that CTW 
scattering strength varied with both bottom topography and 
density stratification creating an extremely complex 
interaction between these processes. The region studied in 
Segment 2 of this analysis appeared to be influenced by 
similar complicated interactions. A model specifically 
designed to account for the highly variable bathymetry of the 
Monterey Canyon would be necessary to produce more meaningful 
results . 

0. THE MONTEREY BAY MEANDER 
1. General Description 

Satellite imagery from the Advanced Very High 
Resolution Radiometer (AVHRR) on the NOAA 11 satellite reveals 
that a distinct meander dominated the surface flow on 01 
September 1991 (Figure 27) . This picture provides the best 
image during the time frame of the energetic feature which 
dominated Segment 3 . The image exhibits characteristics 
similar to those observed by Tracy (1990) . In particular, a 
warm core, anti -cyclonic meander with a radius of 
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Figure 27. AVHRR Sea Surface Temperature Image for 

01 September 1991: The different gray shades 

correspond to different temperatures with white 
representing the coldest water (<12'* C) and 
black the warmest water (17*^ C) . Velocity 
vectors from 400 m at P2 (to the south) and MB2 
(to the north) are shown in white and indicate 
equatorward flow at both moorings. Vectors are 
proportionally scaled for visual effect. 



97 




approximately 65 km can be inferred off Monterey Bay from the 
temperature gradient in the image (Figure 27) . The meander is 
bounded to the south by a pronounced cold filament extending 
some 150 km off Point Sur and to the north by a pocket of cold 
water off Ano Nuevo probably due to intense upwelling often 
observed there. Tracy (1990) described conditions that 
allowed a similar meander to move shoreward and cause radical 
temperature changes inside the Monterey Bay. 

It is important to examine whether features such as 
this impact the currents at 400 - 500 m depth. Ramp et al . 
(1991), Buyer et al. (1991), and Kosro et al. (1991) described 
a similar meander and cold filament off Point Arena in June 
1987. The Point Arena meander displaced isopycnals to depths 
greater than 500 m. Geostrophic velocities calculated by 
Robson (1990) confirmed that surface flow reached depths of 
500 m just outside of Monterey Bay. Additional measurements 
at mooring P3 at other times show strong eddy- like features 
extending to depths > 1000 m during spring 1990 and 1991. The 
current patterns described in Segment 3 also suggested that 
the surface feature extended to the instrument depths at P2 
(500 m) and MB2 (418 m) . 

Poleward currents dominated the flow at both P2 and 
MB2 prior to 07 August (Figures 19 and 21) . It is believed 
that the shoreward movement of the meander caused this 
poleward flow to reverse at P2 first and then MB2 . From 07 to 
24 August, currents were onshore at MB2 and offshore at P2 



98 



suggesting that MB2 was located in the northern boundary of 
the meander while P2 was in the southern boundary. From 14 
August until 03 September, the flow was equatorward at both 
moorings indicating that the surface jet seen in the satellite 
image extended down to the current meter depths at P2 and MB2 . 
The mean current vectors from 01 September at 400 m are 
superimposed on the image, and verify this equatorward flow on 
the edge of the meander. 

Temperature records at P2 and MB2 also showed the 
meander's influence at depth. As the meander moved onshore, 
the temperature at the two instruments increased steadily 
(Figure 28) . After this rise frequent temperature 
fluctuations, perhaps associated with the highly variable 
frontal position, were observed at both moorings (Figure 28) . 
As the meander retreated, both moorings experienced a sharp, 
simultaneous drop in the temperature around 24 August possibly 
due to the cold tongue of water surrounding the meander 
(Figure 27) . 

The Segment 3 vector velocity plot (Figure 19) 
presented another remarkable issue. Shortly after the onset 
of the deep meander near Monterey Bay, steady poleward flow at 
the Farallon moorings D and E, weakened significantly and 
actually reversed to equatorward. Could deep interaction 
nearly 200 km to the south somehow "shut off" the source of 
the poleward flowing California Undercurrent? Although at 
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Figure 28. Temperature Time Series at P2 500 m and MB2 418 

m: Plots show the rise and subsequent drop in 

temperature at the moorings due to across -shore 
excursions of the Monterey Bay Meander. 



first glance there appeared to be some connection, it would be 
imprudent to speculate further without additional evidence. 



E. EXPLAINING THE EOF'S 

1. Mode 1 California Undercurrent 

EOF analysis provided additional evidence to support 
the presence of the CUC. Throughout all three time segments, 
first mode EOF results consistently depicted a poleward flow 
at all locations except MBl, located inshore of the CUC in the 
Monterey Canyon. The utility of this technique relied on the 
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orthogonality principle; i.e., eigenvalues of a covariance 
matrix are orthogonal to each other and allow the flow to be 
separated into distinct patterns. EOF analysis works best 
when the first few modes explain most of the variance. 
Wallace and Dickinson (1972) described this as a necessary 
criteria to establish physical significance. The first mode 
alongshore eigenvalues explained between 64% and 66% of the 
variance with higher modes decreasing significantly 
(Figure 29) . Therefore, the mode 1 result is attributed to 
the dominant poleward flow of the California Undercurrent. 

2 . Segment 2 

The mode 2 EOF structure is most likely due to the 
presence of significant deviations from the mean flow. In 
Segment 2, the mode 2 structure can be attributed to the 
coastal trapping phenomenon believed to dominate the segment. 
Figure 13 (bottom left) shows this idea the best. 
Specifically, the mode 2 structure at P3 opposed that at P2 
and MB2 suggesting that the trapping occurred most efficiently 
shoreward of P3 . 

3 . Segment 3 

In Segment 3, both mode 1 and mode 2 eigenfunctions 
displayed characteristics consistent with the Monterey Meander 
that dominated the flow at P2 and MB2 (Figure 22) . In the 
alongshore component, the mode 2 eigenfunction showed that 
both MB2 and P2 had the same sign which opposed the flow at D 
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Figure 29. Percent Variance Decrease for Alongshore 

Elgenmodes : Shown are the percentages of variance 

explained (ordinate) vs. mode number (abscissa) 
for Segment 1 (top) , Segment 2 (middle) and 
Segment 3 (bottom) . 
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and E. This could be attributed to the strong equatorward 
influence of the meander as described above. In the across - 
shore component, mode 1 showed this structure the best with 
flow at P2 opposing that at MB2 . This can be explained by the 
onshore flow at MB2 opposing the offshore flow at P2 while the 
meander passed through these moorings. 



103 



VI. CONCLUSIONS AND RECOMMENDATIONS 



Three time segments of sub- thermocline (350 - 500 m) 
current meter data along the central California coast were 
examined to determine alongshore low frequency variability on 
time scale of two days to several months. Time and frequency 
domain analyses yielded three basic conclusions: (1) the 
poleward flowing California Undercurrent (CUC) dominated the 
mean flow during all segments; (2) the presence of a coastal - 
trapped significantly influenced the flow between P2 and MB2 
during Segment 2 and (3) the deep flow at P2 and MB2 was 
influenced by a strong equatorward meander of the California 
Current . 

The mean flow in all segments closely resembled that 
expected for the CUC. Segment 1 (December 89 - April 90) most 
nearly followed this mean pattern with few significant 
deviations. The first empirical mode in all segments 
consistently revealed a poleward current structure. Seasonal 
maxima occurred during February and May and showed close 
agreement with other findings in the vicinity of Point Sur. 

Segment 2 (August 90 - March 91) displayed a significant 
deviation from the mean flow in the form of a coastal- trapped 
wave. The observed wavelengths between those for an internal 
Kelvin wave and a continental shelf wave were suggestive of a 
complex interaction between density stratification and the 
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complicated bathymetry of the Monterey Canyon. The wave 
influenced P2 and MB2 much more than P3 , located outside of 
the internal Rossby radius. Though not specifically addressed 
in this research, wave forcing models developed by other 
researchers yielded results not inconsistent with the possible 
forcing mechanism in the Segment 2 wave. Further study of 
coastal - trapping should include an intensive numerical 
analysis that considers both density stratification and 
complex bathymetry. 

Segment 3 (May 91 - November 91) represented another 
deviation from the mean poleward flow. An anti- cyclonic 
meander, clearly visible from satellite imagery, appeared to 
influence the flow down to 500 m. Cross-shore velocity 
component analysis at P2 and MB2 showed a strong negative 
correlation indicating the presence of an anti- cyclonic 
rotating feature. Additionally, temperature time series 
fluctuated in a manner that clearly showed the meander's 
influence. The meander appeared to disrupt the California 
Undercurrent some 200 km to the north at the Farallones 
moorings. An intensive hydrographic study encompassing 
Monterey Bay would provide insight into the surface feature- 
depth interaction and possibly improve the understanding of 
the California Undercurrent forcing. 
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